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ABSTRACT 


This  study  presents  two  methods  of  analysis  for  the 
supersonic  oscillating  cascade  with  subsonic  leading  edge. 

A relatively  straightforward  solution  is  developed  for  the 
slowly  oscillating  finite  and  infinite  flat-plate  cascades 
which  provides  simple  analytical  expressions  for  the  unsteady 
pressure  distributions.  Comparison  with  other  solutions 
is  generally  excellent.  Some  additional  linear  topics 
including  resonance  and  a unique  inflow  condition  are  also 
treated.  In  addition  a nonlinear  method  of  characteristics 
solution  for  finite  cascades  is  described  which  permits 
analysis  of  blade  thickness  effects  on  flutter.  At  this 
time,  only  the'  inlet  and  passage  flow  computations  have 
been  completed  and  are  compared  with  available  experimental 
data. 
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1 . INTRODUCTION 


1.1  General.  Unsteady  flow  in  turbomachines  can  be  broadly 
divided  into  two  categories:  circumferential  asymmetry  in 
the  flow  through  the  turbomachine  and  flow  perturbation  due 
to  the  self-excited  oscillatory  motion  (flutter)  of  the 
internal  blading  [1].  As  shown  in  figure  1.1.1,  taken  from 
[2],  at  least  four  types  of  flutter  may  occur  in  a compressor 
rotor.  The  two  stalled-f lutter  phenomena  are  encountered 
when  compressor  or  fan  blading  is  subjected  to  incidence 
angles  near  or  exceeding  the  stall  angle,  and  so  are  off- 
design  phenomena.  Subsonic  stalled  flutter  has  long  plagued 
the  engine  designer  [3] . Choke  flutter  and  supersonic  unstalled 
flutter  are  low-incidence  instabilities  occurring  when  the 
flow  over  the  blade  is,  respectively,  transonic  along  most 
of  the  chord,  or  supersonic  over  the  outer  span  portions  [3] . 
Supersonic  unstalled  flutter  (hereafter  referred  to  simply 
as  supersonic  flutter) , which  has  appeared  only  recently  as 
a significant  problem  in  turbomachinery,  is  the  subject  of 
this  investigation.  Primarily  a concern  in  compressors  and 
fans  [1] , it  is  an  instability  born  of  the  recent  technologi- 
cal advances  aimed  at  increasing  the  thrust-to-weight  ratios 
of  turbomachines:  increased  tip  speeds,  composite  blades, 
lighter  discs  and  blades,  and  the  elimination  of  part-span 
shrouds  [2,3].  Some  insight  into  the  nature  of  supersonic 
flutter  is  gained  by  a consideration  of  its  characteristics: 


Figure  1.1.1  Compressor  map  showing  boundaries  for  four 

types  of  flutter. 

1.  Its  flutter  boundary  imposes  a high-speed  limit  on 
compressor  operation  [2] . 

2.  It  is  a low- incidence,  attached-flow  instability  that 
can  occur  at  the  design  condition  (figure  1.1.1). 

3.  Flutter  amplitude  decreases  with  increasing  back 
pressure  on  the  rotor  [2].  (For  this  reason  the  present 
analysis  considers  the  most  critical  case  of  zero  pressure 
rise  across  the  rotor.) 

4.  At  the  flutter  condition  all  the  blades  are  observed 


to  flutter  at  the  same  frequency  with  a constant  interblade 
phase  angle  [2]. 


5.  Both  single-degree-of-freedom  torsional  flutter  (pitch) 


and  two-degree-of-f reedom  torsional-flexural  flutter  (pitch 
and  plunge)  can  occur  [3,4]^ 

6.  The  flutter  boundaries  are  significantly  affected  by 
airfoil  shape  [4,5]. 

7.  The  primary  flutter  mechanism  is  an  inviscid  phenomenon. 
If  the  unsteady  aerodynamic  force  lags  the  blade  motion 

by  a phase  angle  which  is  sufficient  to  ensure  positive 
work  over  the  entire  cycle,  then  flutter  can  occur  [3,4]. 

It  is  clear  that  reliable  techniques  for  predicting  the 
unsteady  pressure  distribution  on  a fluttering  blade  are 
required  if  one  hopes  to  predict  and  avoid  the  onset  of 
flutter.  With  this  end  in  mind  the  cascade  has  been  exten- 
sively employed  as  a theoretical  abstraction  of  a compressor 
rotor.  The  unsteady  flow  is  considered  in  an  annulus  of 
differential  radial  height  within  which  equispaced  airfoils 
are  radially  disposed.  Unwrapping  such  an  annulus  gives 
rise  to  the  two-dimensional  approximation  of  plane  flow 
through  a cascade  of  airfoils  [1] . The  cascade  is  termed 
linear  if  the  airfoils  are  flat  plates  of  zero  thickness. 

As  seen  in  figure  1.1.2,  a rotor  does  not  require  a super- 
sonic axial  flow  to  have  a supersonic  velocity  component 
relative  to,  and  parallel  to,  the  blade  axis.  Indeed  this 


^R.A.  Arnoldi  ejt  in  "Supersonic  Chordwise  Bending 

Flutter  in  Cascades,"  Pratt  & Whitney  Report  No.  PWA-5271 
of  31  May  1975  show  that  chordwise  bending  can  also  occur. 


f 


Figure  1,1.2  Supersonic  cascade  classification. 


distinction  provides  a convenient  method  of  classifying 
supersonic  cascades.  If  the  axial  velocity  is  subsonic, 
the  cascade  is  said  to  have  a subsonic  leading  edge;  in  the 
other  case  it  has  a supersonic  leading  edge  [6],  The  physi- 
cal significance  of  this  distinction  is  clear  from  figure 
1.1.2.  In  the  supersonic-leading-edge  cascade  each  blade 
influences  only  its  immediate  neighbors.  When  the  cascade 
has  a subsonic  leading  edge,  however,  disturbances  generated 
by  the  upstream  portions  of  each  blade  affect  the  flow  around 
each  subsequent  blade  in  the  cascade;  furthermore  the  flow- 
field  at  the  underside  of  each  blade  is  influenced  by  the 
wake  disturbances  of  all  preceding  blades.  Since  the  supersonic- 
leading-edge  cascade  has  been  thoroughly  investigated  [7,8,9, 

10],  the  more  difficult  subsonic-leading-edge  case,  which  is 
also  the  more  interesting  in  practice,  is  considered  in  the 
present  study.  Another  cascade  classification  is  frequently 
employed;  if  each  blade  is  assumed  to  be  preceded  and  followed 
by  infinitely  many  others,  the  cascade  is  said  to  be  infinite; 
if  there  is  an  identifiable  first  blade  — a necessary  pre- 
requisite for  many  numerical  solutions  — the  cascade  is 
termed  finite. 

A considerable  amount  of  effort  has  been  expended  in 
attempting  to  predict  — numerically  and  analytically  — 
the  unsteady  pressure  distribution  over  a blade  in  an  oscilla- 
ting supersonic  cascade  with  subsonic  leading  edge.  A 
brief  discussion  of  some  of  the  more  noteworthy  advances 


follows . 


1 . 2 The  Linear  Subsonic-Leading-Edge  Cascade.  Verdon 


and  McCune  [11]  formulated  a linear  boundary  value  problem 
for  the  subsonic-leading-edge  configuration,  solving  it 
for  velocity  potential.  This  permitted  calculation  of  the 
pressure  along  the  entire  upper  surface  of  an  arbitrary  blade 
and  along  the  lower  surface  from  the  leading  edge  to  the 
point  of  impingement  of  the  trailing-edge  Mach  wave  from 
the  preceding  blade.  In  a more  recent  paper  Verdon  [12] 
formulated  a second  boundary  value  problem,  in  which  pressure 
is  the  basic  dependent  variable,  which  permitted  determination 
of  the  pressure  along  the  entire  lower  surface  as  well. 

Nagashima  and  Whitehead  [13]  and  Goldstein  [14]  obtained 
solutions  by  replacing  the  cascade  blades  with  doublet  dis- 
tributions. Kurosaka  [6]  introduced  what  he  terms  the 
passage  approach  in  applying  Laplace  transform  techniques 
to  a slowly-oscillating  cascade,  obtaining  closed-form 
analytical  expressions  for  the  pressure  distributions  of 
certain  configurations.  He  recently  [15]  extended  his  solu- 
tion to  include  arbitrary  frequencies  of  oscillation.  In 
assuming  that  a space  periodicity  of  the  flow  exists  from 
one  blade  passage  to  the  next,  Kurosaka  claims  that  he 
avoids  complications  arising  from  the  breakdown  of  linear 
acoustic  theory  in  the  far  field  [4] . Platzer  and  Brix 
[16]  and  Platzer  and  Bell  [17]  employed  the  method  of  char- 
acteristics to  solve  directly  for  the  perturbation  veloci- 
ties and  pressures  everywhere  in  the  flowfield.  Their 
computer  program  has  been  thoroughly  checked  against  existing 

^ 


Figure  1.2.1  Limits  of  applicability  of  Chadwick's  solution. 


solutions  [11,12,13]  and  is  extensively  employed  in  the  ! 

i 

present  study  as  a method  of  comparison  for  both  the  analy- 
tical and  numerical  solutions  presented.  Chadwick  [5]  used 

i 

Sauer's  solution  [18]  of  the  unsteady  potential  equation, 
accurate  to  a first  order  in  frequency  of  oscillation,  to 
determine  the  pressure  distribution  on  a given  blade  from 
the  leading  edge  to  the  point  of  the  first  reflection  (point 
A in  figure  1.2.1)  of  disturbances  from  the  preceding  blade. 

He  obtained  exact  analytical  agreement  with  Kurosaka  [6]  over 
that  portion  of  the  blade  for  the  single-degree-of-freedom 
pitching  cascade.  Chadwick's  technique  has  the  decided 
advantage  of  being  much  simpler  than  that  of  Kurosaka,  easily 




yielding  closed-form  expressions  for  any  cascade,  and 
raising  the  question  as  to  whether  the  method  could  be 
extended  to  the  entire  blade  for  the  cases  of  both  pitch 
and  plunge.  Finally,  Adamczyk  [19]  , in  a yet-unpublished 
work,  employs  the  Wiener-Hopf  method  to  obtain  the  pressure 
distribution  for  the  linear  cascade  with  subsonic  leading 
edge. 


1.3  The  Nonlinear  Subsonic-Leading-Edge  Cascade.  Attempts 
to  deal  analytically  or  numerically  with  an  oscillating 
cascade  of  blades  having  thickness  are  not  nearly  so  numerous 
as  their  linear  counterparts.  Carrier  [20]  provided  an  exact 
theory  for  a single  oscillating  wedge.  VanDyke  [21,22]. 
also  treated  the  isolated  wedge  as -well  as  an  oscillating 
airfoil  of  arbitrary  profile.  Teipel  [23,24]  employed  an 
ingenious  method-of-characteristics  approach  to  numerically 
obtain  the  unsteady  pressures  on  a single  oscillating  profile 
of  arbitrary  shape,  obtaining  good  agreement  with  the  single- 
blade solutions  of  Carrier  and  VanDyke.  Teipel 's  method 
assumes  potential  flow  with  negligible  entropy  changes 
across  leading-edge  and  trailing-edge  shocks.  Kurosaka  and 
Edelfelt  [4]  have  obtained  numerical  results  for  an  isolated, 
oscillating,  parabolic-arc  airfoil  and  are  presently  working 
to  extend  their  solution  to  a cascade  of  such  airfoils. 
Chadwick  [5]  employed  Teipel 's  method  to  obtain  the  unsteady 
pressure  distribution  over  the  preinterference  zone  of  an 
arbitrary  blade  in  a cascade  of  wedges.  His  results  indicate 


that  nonlinear  thickness  effects  significantly  alter  the 
unsteady  pressure  profile  of  an  arbitrary  blade. 

1.4  The  Present  Investigation.  The  present  work  has  two 
major  parts:  The  first  is  analytical  and  treats  the  linear 
cascade  while  the  second  is  numerical,  dealing  with  a non- 
linear configuration. 

In  the  linear  effort  a closed-form  analytical  solution 
for  the  pressure-difference  distribution  along  a blade  in 
a slowly-oscillating  cascade  is  obtained.  The  solution 
technique  is  that  developed  by  Chadwick  [5]  in  the  pre- 
interference zone  (figure  1.2.1)  of  a pitching  cascade,  in 
which  Sauer's  solution  of  the  unsteady  linearized  potential 
equation  is  employed.  In  the  present  work  Chadwick's  solu- 
tion is  extended  through  all  reflection  zones  and  into  the 
blade's  wake.  In  addition,  the  solution  is  obtained  for 
two  degrees  of  freedom  (pitch  and  plunge) . Two  standard 
configurations  — Verdon's  cascades  A and  B [12]  — are  con- 
sidered, one  from  a finite-cascade  approach  and  the  other 
as  an  infinite  cascade  in  which  Kurosaka's  passage  approach 
is  employed.  Extensive  comparisons  are  made  with  existing 
solutions  [6,12,19,17];  and  several  special  topics,  including 
cascade  resonance  and  a unique  inflow  condition,  are  briefly 
treated. 

The  nonlinear  effort  is  an  evaluation  of  the  feasibility 
of  extending  Teipel's  unsteady  method  of  characteristics 
to  a cascade  of  airfoils  with  arbitrary  profiles.  Teipel 


[23,24]  represents  the  nonlinear  flowfield  by  a harmonic, 
unsteady,  velocity  potential  superimposed  on  a steady  one, 
and  discontinuities  in  the  flowfield  variables  across  the 
shockwaves  emanating  from  the  airfoil  are  determined  from 
the  isentropic  Rankine-Hugoniot  relations.  These  concepts, 
coupled  with  the  characteristic  grid  and  compatibility  relations 
derived  by  Teipel,  permit  determination  of  the  steady  and 
unsteady  preinterference  pressures  acting  on  the  second 
blade  of  a cascade  of  parabolic-arc  airfoils.  The  numerical 
results  are  compared  with  experimental  data  from  tests  con- 
ducted with  a pitching  cascade  by  Fleeter  and  Riffel  [25.] 
and  with  linear  theory.  An  improved  method  for  determining 
the  jump  conditions  across  the  bow  shock  is  presented  and 
the  computatio-nal  procedure  is  carefully  outlined  to  facili- 
tate extension  of  the  technique  to  the  entire  flowfield  of 
a completely  general  cascade. 

1.5  Some  Comments  on  Notation.  In  the  linear  portion  of 
this  investigation  it  is  the  unsteady,  dimensionless, 
pressure-difference  distribution  along  a particular  cascade 
blade  that  is  sought.  This  pressure  difference  is  a complex 
quantity  given  the  symbol  and  represents  an  amplitude 

function  of  distance  along  the  blade.  The  harmonic  oscilla- 
tion  factor  e is  always  implicitly  understood.  In 

addition  it  is  desirable  to  obtain  an  expression  for 


ACp  at  the  instant  that  the  blade  in  question  achieves 


maximiim  up  displacement,  referred  to  as  the  initial  position. 

This  expression  is  obtained  by  multiplying  through  by  the 

factor  e and  so  it  is  indicated  by  the  sumbol  AC  e . 

P 

In  the  nonlinear  portion  of  this  study  only  the  local 


pressures  — Cp  upper  and  lower  — are  sought.  They  are 
also  complex  quantities  and  it  is  convenient  to  express 
them  as  a magnitude  and  phase  angle  relative  to  the  motion 
of  the  blade  under  study;  the  symbols 
respectively,  are  employed. 
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2 . THE  LINEAR  LOW-FREQUENCY  CASCADE 


2.1  General.  A cascade  of  flat  plates  with  chord  c 
oscillating  harmonically  with  frequency  in  a steady, 
uniform  flow  U is  said  to  oscillate  at  a dimensionless 
reduced  frequency  of  oscillation  k defined  by 

(2.1.1) 

2 

If  k is  small  when  compared  with  k , then  the  cascade  is 
termed  a low-frequency  or  slowly-oscillating  cascade.  Although 
low-frequency  flutter  is  not  common,  it  is  worthy  of  investi- 
gation for  se.veral  reasons,  the  major  one  being  that  it  is 
the  only  case  for  which  one  can  hope  to  attain  a relatively 
simple,  closed-form  expression  for  the  unsteady  pressures 
on  a given  blade.  Such  expressions  provide  insight  into  the 
flutter  phenomenon,  as  evidenced  by  the  identification  of 
a unique  inflow  condition  in  a subsequent  chapter  of  this 
study.  Also,  the  low-frequency  cascade  provides  an  impo -tant 
limiting  case  for  any  general-frequency  solution,  numerical 
or  analytical. 

2 . 2 Problem  Formulation. 

2.2.1  The  Linearized  Potential  Equation.  The  flow 
varicdales  in  the  cascade  flowfield  are  assumed  to  be 
composed  of  a freestream  component  and  a much  smaller 


perturbation  component 


(2.2.1) 


The  cascade  is  assumed  immersed  in  a perfect  gas,  each 
particle  of  which  maintains  constant  entropy  as  it  moves 
through  the  flowfield 


D-t 


s=0 


(2.2.2) 


In  addition  the  flow  is  assumed  inviscid  and  irrotational, 
guaranteeing  the  existence  of  a velocity  potential 

(2.2.3) 


The  continuity  and  momentum  equations  for  the  flow  are, 
respectively 


V*  V = o 


pv 

D-fc. 


(2.2.4) 

(2.2.5) 


introducing  the  velocity  potential  and  integrating  (2.2.5) 
produces,  after  linearization 


at 


Px 


+ 


o 


(2.2.6) 


Using  the  equation  for  the  freestream  speed  of  sound 


(2.2.7) 


the  continuity  equation  can  be  line.arized  to  give 


aV  \ 

5^*  /o.  oi 


(2.2.8) 
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2,2.2  The  Linearized  Pressure  Coefficient.  A generalized 
form  of  Bernoulli's  ecuation  follows  directly  from  the 
momentum  equation  (2.2,5) 
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:s  CONST 


(2.2.13) 


- 2 


Introducing  the  perturbation  quantities  in  the  form  of 
(2.2.1)  and  ignoring  the  products  of  perturbation 
quantities  one  obtains 


If  +xj|f  + ^=o 

00 


(2.2.14) 


where  the  prime  notation  has  again  been  dropped.  Thus,  the 
dimensionless  pressure  coefficient  is  given  by 


r - 'P-'P* 


2.2.3  The  Linearized  Boundary  Conditions.  The  surface 
of  a body  placed  in  the  flow  can  be  described  as  a function 
of  space  and  time 


o 


(2.2.16) 


where  r is  the  position  vector  of  points  on  the  surface. 

The  condition  that  the  flow  at  the  body  surface  be  everywhere 
tangent  to  the  surface  can  be  expressed  by 


(2.2.17) 
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Figure  2,2.1  Two-degree-of-freedom  oscillation. 

at  the  surface.  If  the  body  is  a flat  plate  oscillating 
harmonically  in  plunge  and/or  pitch,  as  shown  in  figure  2.2.1, 
then  the  describing  functions  are  respectively 

(2.2.18) 

(2.2.19) 

where  h^  and  are  the  maximum  displacements  and  b 

is  the  pitch  axis  location  on  the  x axis . Applying  the 
generalized  flow  tangency  condition  (2.2.17)  to  (2.2.18) 
and  (2.2.19)  yields  the  boundary  conditions  for  plunge  and 
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pitch  respecitvely 


= -ho  a.-t  ^ = o (2.2.20) 

*1^  = c-t^r  O (2.2.21) 

2.2.4  Nondimensionalization . It  is  convenient  to  intro- 
duce nondimensional  variables  into  the  governing  equations 
to  simplify  and  generalize  them.  The  distance  variables 
(x,y,b)  are  normalized  by  chordlength  c , velocity  potential 
by  the  product  Uc  and  time  by  the  quantity  c/U  . It  is 
also  convenient  to  set 

h,=  c 1.  RADIAN  (2.2.22) 

and  to  separate  the  velocity  potential  into  a complex 
amplitude  and  phase  angle 

(2.2.23) 

Thus  the  potential  equation  (2.2.11),  pressure  coefficient 
(2.2.15)  and  boundary  conditions  (2.2.20,21)  become 

(2.2.24) 

— (2.2.25) 





r 


Co=  - 
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where  ^ 


(2.2.28) 

and  Cp  and  Jjl  are  complex  amplitude  functions  of  x and 

ilcti 

y only  with  the  factor  e omitted. 


2.2.5  Cascade  Geometry.  The  physical  parameters  chosen 
to  describe  the  cascade  are  the  vertical  distance  between 
the  blades’  d and  the  horizontal  distance  between  the 
leading-edge  Mach  waves  m . Both  parameters  are  normalized 
by  chordlength.  The  cascade  can  also  be  described  by 
solidity  fi  and  stagger  angle  0 , which  are  shown  in 
figure  2.2.2.  These  parameters  are  related  as  follows 


(2.2.29) 

(2.2.30) 


where  y8  is  given  by  (2.2.28).  The  linear  cascades  employed 
in  this  study  are  Verdon's  cascades  A and  B [12],  which 
have  become  standards  in  the  field  and  are  depicted  in 
figure  2.2.2. 


Si 


2.3  The  Isolated  Airfoil.  Consider  a single  flat  plate 
oscillating  with  two  degrees  of  freedom  as  shown  in  figure 
2.2.1,  and  define 

(2.3.1) 

Then  s is  constant  along  any  left-running  (upper  surface) 
Mach  wave,  while  s is  constant  along  right-running  (lower 
surface)  waves.  Sauer  [18]  showed  that  ti\e  potential 
equation  (2.2.24)  has  the  following  left-running  and 
right-running  solutions,  respectively,  for  a slowly- 
oscillating  airfoil 

(2.3.2) 

(2.3.3) 

where  the  g and  h functions  are  determined  from  the 
boundary  conditions  (2.2.26,27).  Requiring  that  the 
velocity  potential  be  continuous  across  the  leading-edge 
Mach  waves,  one  has 

^^o(o^=  V»„(o')  = sV,^to^  so  (2.3.4) 


Enforcing  flow  tangency  at  the  upper  surface  of  the  airfoil 
yields,  for  plunge  and  pitch  respectively 


Figure  2.2.2  Cascade  geometry. 


PL:  = 0 


>,„CSJ  = ~ 


(2.3.5) 


PI: 


ViulS^s-^j(  ^-*-byj*s)  (2.3.6) 


And  so  on  the  upper  surface 


PL: 


PI: 


(2.3.7) 


8) 


Thus  the  upper-surface  pressure  coefficient  follows  from 
(2.2.25) 


PL: 


PI: 


.u  _ a.1^ 

r--7' 


(2.3.9) 


Cls 


(2.3.10) 


where  the  subscript  indicates  the  single  airfoil.  Similarly, 
for  the  lower  surface 


(2.3.11) 


2.4  Continuity  of  the  Velocity  Potential  Across  a Mach  Wave. 
Consider  two  reflection  zones  separated  by  the  Mach  wave 
s = s^  , as  shown  in  figure  2.4.1,  and  let  the  variables  z 
and  z respectively  measure  distance  tangent  and  normal  to 
the  Mach  wave.  The  velocity  potential  in  zone  II  is  simply 
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Figure  2.4.1  Variation  of  the  velocity  potential  along 
the  boundaries  of  reflection  zones. 


that  in  zone  I plus  an  additional  disturbance 


(2.4.1) 


Since  the  tangential  velocity  component  remains  unchanged 
across  the  wave,  (2.4.1)  implies  that  is  a function  of 

z alone,  along  the  Mach  wave 


•I'OX  &e  (2.4.2) 


But  since  z is  itself  constant  along  the  Mach  wave,  then 
must  have  the  form 


CONST.  4or 


(2.4.3) 


J 


That  the  constant  is  zero  follows  from  the  physical  con- 
straint that  the  velocities  remain  finite  in  the  wave 
itself,  and  so 

IS*.-)  (2.4.4) 

which  is  in  agreement  with  the  conclusion  reached  by 
Landahl  [26]. 

2 . 5 The  Finite  Two-Degree-of-Freedom  Cascade.  As  an 
elementary  means  of  introducing  the  method  of  analysis, 
the  pressure-difference  distribution  is  found  for  the  second 
blade  of  a finite  version  of  cascade  A.  This  analysis  also 
affords  the  opportunity  to  compare  the  finite-cascade  results 
obtained  with  numerical  finite-cascade  data  [17]  and  analy- 
tical results  [6]  of  existing  solutions.  It  is  emphasized 
that  the  pitch  results  through  equation  (2.5.27)  were  first 
obtained  by  Chadwick  [5] . The  reflection  zones  referred  to 
are  pictured  in  figure  2.5.1,  and  the  pressure  results  are 
given  in  terms  of  the  isolated-blade  pressure  coefficients 
(2.3.9,10).  Before  proceeding  with  the  analysis,  some 
additional  variadales  are  introduced  for  blade  two 

(2.5.1) 

which  are  related  to  those  of  blade  one  by  the  translation 
relations 
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Cascade  A (not  to  scale) 
Mach  nvunber:  1.345  / 


blade  3 


blade  2 ^ / 


-IIIA 


IV  \ V \ VI 


bl?de  1 


Figure  2.5.1  The  finite-cascade  reflection  zones. 


(2.5.2) 


S=C'-^m 


Zysa  (2.5.3) 


2.5.1  Zones  I and  lA.  In  zones  I and  lA  Sauer's  single- 
blade solution  applies 


S6f  (2.5.4) 

(2.5.5) 


and  in  zone  lA  the  boundary  conditions  take  the  form 
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PI:  -|^  = at.  1*0  (2.5.7) 


where  yw.  is  the  interblade  phase  angle  by  which  the  second 
blade  leads  the  first.  Applying  the  boundary  condition  at 
the  upper  surface  of  each  blade  and  using  that 

= W.to')  =0  (2.5.8) 

gives,  for  plunge  and  pitch 


(2.5.9) 

= ^ B (2.5.10) 

-^G  (2.5.11) 

F 

(2.5.12) 

-4-  bys’-o'i  ev*  G M*(n  4-  j> 

where  A,  B,  C and  D are  constants  of  integration.  By 
enforcing  continuity  of  the  velocity  potential  across  the 
leading-edge  Mach  wave  of  blade  two,  the  constants  are 
evaluated  and  the  velocity  potentials  determined  in  terms 
r of  the  isolated  blade  potential  ?ir(2-3.7.8) 

L 


PL: 


Pi! 


s O 


PL: 

r 


PI: 


h^CCT) 


= ->[c¥ 


where  (2  is  a constant  that  need  not  be  evaluated  because 

2 


it  becomes  the  coefficient  of  a term  in  k when  the 

1 


pressure  coefficients  are  formed.  Thus,  on  the  upper 
surface  of  the  second  blade,  in  terms  of  the  isolated-blade 
pressure  coefficient  Cp,  (2.3.9,10) 


PL; 


= Cp- 


(2.5.15) 


PI; 


P 


(2.5.16) 


2.5.2  Zones  IIB  and  VIB.  In  zone  IIB  the  velocity 
potential  is  that  of  the  isolated  blade  (2.3.3) 


{2.S.11) 


while  in  VIB  the  potential  acquires  an  additional 


The  symbol  0 is  used  throughout  section  2 to  indicate 
a constant  which  need  not  be  determined.  *• 
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right-running  disturbance 


<^1®*  (2.5.I8) 


2.5.3  Zone  II.  In  zone  II  a right-running  reflection 
occurs  and,  using  the  translation  relations  (2. 5. 2, 3)  to 
shift  to  the  blade-two  coordinate  system,  one  has 


PL;  9^?  = (2.5.19: 


(2.5.20) 


where 


=s  V>j(0^  = O 


(2.5.21) 


Applying  the  flow  tangency  condition  at  the  lower  surface 
of  blade  two  and  using  the  Sauer  functions  to  construct  the 
velocity  potentials  yields 


(2.5.22) 

0®- A'-*- 

•V  (2.5.23) 


And  the  pressure  coefficients  follow  from  the  potentials 


PL:  = Cr,(2-e^^)  (2.5.24) 

PI:  Cpf- =:  Cp,  (2.5.25) 

Thus  the  pressure  difference  can  be  computed  through  zone 
II,  and  by  multiplying  by  the  factor  e rhe  pressure- 

difference  amplitude  is  obtained  at  the  moment  the  second 
blade  reaches  maximum  displacement 

PL:  i^Cpe''-^  0±  ^ < \-m  (2.5.26) 

ACpe-‘>»*  = 

PI:  (2.5.27) 

o i ^ <.  \-m 

2.5.4  Zone  III.  In  zone  III  a left-running  wake  distur 
bance  is  experienced,  and  after  applying  the  translation 
relations 

PL:  9^7**^  - (2.5.28) 

- S (Vwi  M*)  [-yflhM*-*-  - e'/*)]  (2.5.29) 

-*•^1 


where 


^jCO  = h,Cl)  = O 


(2.5.30) 


Since  zones  III  and  VIB  are  separated  by  the  wake  of  the 
blade  rather  than  the  blade  itself,  the  wake  boundary  con- 
ditions are  used  between  the  two  zones.  The  first  of  these 
demands  pressure  equilibrium  across  the  wake 


_ -SIB 


O-t. 


so 


(2.5.31) 


which  leads  to  the  equations 


= ST  (Cp,) 


(2.5.32) 


V>3  + 1 = FR  (.  C.l'i  (2.5.33) 


where  the  right-hand  sides  of  (2.5.32.-33)  refer  to  the 


steady  and  frequency  terms  respectively  of  C 


u 

pi 


(2.3.9,10) 


The  second  wake  condition  requires  continuity  of  the  normal 
velocity  component  through  the  wake 


(2.5.34) 


which  leads  to  the  equations 


(2.5.35) 
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\x\  (S>  h. 


(2.5.36) 


Equations  (2.5.32,33,35,36)  are  solved  for  the  zone-III 
disturbances 


^j(S')r=0  K;(0)  = --^ 


(2.5.37) 


V{3(s>=^j(sy‘^.b^»-)  (2.5.38) 


2.5.5  Zone  V.  In  zone  V an  unknown  disturbance  is  added 


(2.5.39) 


which  is  determined  by  enforcing  flow  tangency  at  the  lower 
surface  of  blade  two  in  zone  V 


Ws(^>  - 


(2.5.40) 


^yCcr')  = — I -vy>0 


(2.5.41) 


which  leads  directly  to  the  pressure  difference  for  the 
last  portion  of  the  blade 


« 


w 


PL: 


PI; 


(2.5.42) 


(2.5.43) 


l-vn  <.^  t I 


2.5.6  Results  and  Conclusions.  The  finite  results  for 
the  second  blade  are  compared  with  four  existing  solutions 
in  figures  2.5.2  and  2.5.3  for  three  interblade  phase  angles: 
Kurosaka's  solution  [6]  is  a closed-form,  analytical, 
infinite-cascade  solution;  Platzer's  [16,17]  is  a numerical, 
method-of-characteristics  solution  for  the  finite  cascade; 
and  Adamczyk's  [19]  and  Verdon's  [12]  solutions  are  numeri- 
cal, infinite-cascade  results.  It  is  interesting  to  note 
that  the  finite  results  agree  exactly  with  Kurosaka's 
infinite  solution  over  the  entire  blade  for  the  plunge  case 
and  through  zone  II  for  the  pitch  case.  Over  the  zone-V 
portion  of  the  pitching  blade  — the  portion  exposed  to  the 
wake  disturbances  from  blade  one  — the  two  solutions  agree 
numerically,  as  seen  in  figure  2.5.3,  but  are  not  analytically 
the  same.  In  fact,  Kurosaka's  solution  is  singular  for 
yw  = 0 (in-phase  oscillation  of  the  blades)  over  the  last 
portion  of  the  blade,  while  equation  (2.5.43)  behaves  well 
for  that  case.  All  five  solutions  were  found  to  be  in 
excellent  agreement  for  interblade  phase  angles  larger  than 
15®,  but  for  small  angles  considerable  disagreement  is 
evident,  as  seen  in  figure  2.5.2. 
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Second  blade  of  Cascade  A 


M:  1.345 
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Pigiire  2.5.2  Comparison  of  the  finite-cascade 
solution  with  existing  solutions, 
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Figure  2.5.3  Comparison  of  the  finite-cascade 
solutions  with  existing  solutions, 
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Chadwick  [5]  compared  his  finite  and  infinite  results 
for  the  pitch  case  over  the  in-flow  portion  of  the  blade  and 
noted  that  although  the  finite  and  infinite  pressure- 
difference  (ACp)  distributions  are  the  same,  the  local 
pressure  distributions  (C^,  C^)  may  in  fact  differ.  The 
author  found  this  phenomenon  occuring  over  the  remaining 
portion  of  the  pitching  blade  also.  Figure  2.5.4  shows  that 
the  local  pressure  distributions  for  the  present  finite 
solution  and  Verdon's  infinite  solution  do  not  agree  as 
well  as  the  pressure-difference  distributions.  Some  con- 
clusions can  be  drawn  from  the  analysis  at  this  point: 

1.  For  low  frequencies  of  oscillation  the  finite  cascade 
rapidly  approaches  the  infinite-cascade  solution,  achieving 
good  agreement  with  the  pressure-difference  distribution 
of  the  latter  at  the  second  blade. 

2.  For  small  interblade  phase  angles,  the  present  solution 
agrees  with  Kurosaka's  Laplace  transform  solution  for  the 
infinite  cascade  but  considerable  differences  with  the  other 
published  solutions  are  noted. 

3.  Although  the  present  finite-cascade  solution  and 
Verdon's  infinite-cascade  solution  for  the  pressure-difference 
distributions  are  in  agreement  for  large  interblade  phase 
angles,  the  local  pressure  distributions  resulting  from  the 
two  solutions  are  seen  to  disagree. 

2 . 6 The  Infinite  Two-Degree-of-Freedom  Cascade.  In  order 
to  demonstrate  the  applicability  of  the  above  analysis  to 

i 


Figure  2.5.4  Comparison,  of  the  finite  and  infinite 
.*  local  pressure  distributions. 

the  infinite  cascade,  it  is  applied  to  an  infinite  version 
of  Cascade  B.  Since  the  advantage  of  an  identifiable  first 
blade  is  lost,  a condition  specifying  the  space  periodicity 
of  the  flowfields  between  the  blades  is  derived.  The  pitch 
results  through  equation  (2.6.17)  were  first  obtained  by 
Chadwick  [5] , and  the  reflection  zones  referred  to  in  the 
solution  are  depicted  in  figure  2.6.1. 


I 

i 

< 

I 

I 
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Figure  2.6.1  The’ infinite-cascade  reflection  zones. 


2.6.1  The  Flow  Periodicity  Condition.  The  periodicity 
of  the  flow  from  blade  passage  to  blade  passage  was  first 
identified  by  Kurosaka  [6]  and  is  the  heart  of  his  passage 
approach  to  the  infinite  cascade.  An  alternate  derivation 
follows. 

Consider  the  corresponding  zones  I and  lA,  and  assume 
that  the  potentials  are  the  same  there  except  for  a complex 
factor  K 

K (2.6.1) 


Since  each  potential  must  satisfy  the  boundary  condition  of 
its  corresponding  blade 


-Lte  <»> 


«x-t  so 


(2.6.2) 

(2.6.3) 


it  is  clear  that 


K.  = e 


(2.6.4) 


and  so 

(2.6.5) 

The  same  result  follows  for*  any  pair  of  corresponding  zones 
on  adjacent  blades  — for  pitch  and  plunge  — and  is  termed 
the  flow  periodicity  condition. 

2.6.2.  Zones  I and  lA.  The  analysis  in  zones  I and  lA 
proceeds  much  as  in  the  finite  case  (section  2.5.1).  The 
velocity  potentials  have  the  form 

(2.6.6) 

(2.6.7) 

and  the  boundary  conditions  yield 


PL: 


^.(S^rrA  ^,a(^)=.B 


(2.6.8) 


PI;  ^ (2.6.11) 

h,^tar)=  - ^ |^c~  -*•  W^>tr)  e'*--^  -y  m^1>o^  ^ 


The  constants  A through  D are  evaluated  using,  respectively, 
continuity  of  the  velocity  potential  across  the  leading-edge 
Mach  wave  of  the  (n+l)st  blade  and  the  flow  periodicity 
condition 


(2.6.12) 


n+-i 


(2.6.13) 


The  pressure  coefficients,  in  terms  of  those  of  the  isolated 
blade  (2.3.9,10),  follow 


PL; 


(2.6.14) 
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(2.6.15) 
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2.6.3  Zones  II  through  IV.  Finding  the  additional 
reflections  through  zone  IV  is  accomplished  by  routine  use 
of  the  flow  tangency  condition  (2.2.26,27),  and  so  only  the 
pressure  results  for  the  (n+l)st  blade  are  given  below 


^ (2.6.16) 


o ^ ^ <2.^J 


(2.6.17) 


< \-m  (2.6.18) 


' / 


(2.6.19) 


2.6.4  Zones  V,  VI,  and  VII.  In  the  first  three  zones 
above  the  wake  of  the  nth  blade 

[hsU>-  (2.6.20) 

^ (2.6.21) 

9^11^  * (2.6.22) 


where  the  g^  and  h^  Sauer  functions  are  known  from  the 
analysis  in  zone  IV  and 


■ LU  = Vis  CO  :s  O 


(2.6.23) 


>»,C»-Vh)=rO 


(2.6.24) 


2.6.5  Zones  VA  and  VIII.  In  the  zones  straddling  the 
wake  of  the  (n+l)st  blade,  using  flow  periodicity  and 
results  from  zones  IV  and  V 

9?^  (2.6.25) 


where 


5"n-v« 


~ +j^sCcr)e'^>* 

tor) &'">*+  h*C^)  (2.6.26) 

- [<^5(<r)e->-^,Co‘)-^gti')]  - ^4^5(0-)^'^" 


-K IR  [hsr^  (tr)  - 7 (o-)]  = 


(2.6.27) 


CO  = haCl)  =r  o 


(2.6.28) 


The  wake  boundary  conditions  — pressure  equilibrium  and 
normal  velocity  continuity  — are  enforced  between  zones 


VA  and  VIII 


at  " d t 


ok,t.  t=o 


(2.6.29) 


(2.6.30) 


giving  rise  to  two  equations  in  three  unknown  Sauer 


functions 


= ST  (2.6.31) 


(2.6.32) 


+ ^ [?S  I'-  -^7  ->]  = PR  [- 


(2.6.35) 


+ (2.6.34) 


where  the  right-hand  sides  of  (2.6.31,33)  indicate  the 
steady  and  frequency  terms,  respectively,  of  the  enclosed 
quantity.  The  third  equation  in  the  set  comes  from  applying 
the  boundary  condition  to  n+1  lower  surface  of 


the  (n+1) St  blade 


=o 


(2.6.35) 


(2.6.36) 


Solving  equations  (2.6.31)  through  (2.6.36)  for  the  Sauer 
functions  permits  determination  of  the  pressure  difference 
over  the  last  portion  of  the  blade 


PL: 


PI: 


A.Cpe’^>*  = -2.CJJ  »-w<^  < 


^ Cp  e.'>  :s  -v 


^\M.  y»\ 


j-vn<^  <Z.^A  -tvw 


(2.6.37) 

(2.6.38) 
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(2.6.39) 


= -aCp,(a-e.''-'*')  - ^a  (^<1m 
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(2.6.40) 


2.6.6  Results  and  Conclusions.  Figure  2.6.2  shows  that , 
as  with  the  finite-cascade  solution,  there  is  considerable 
disagreement  among  the  several  existing  solutions  for  small 
interblade  phase  angles , but  even  for  those  angles  the 
present  solution  and  Kurosaka's  are  in  agreement.  The 
pressure-difference  expressions  derived  above  for  the  infinite 
cascade  agree  exactly  with  those  of  Kurosaka  [6]  over  the 
entire  plunging  blade  and  up  to  for  the  pitching 

one.  Over  the  last  portion  of  the  pitching  blade  the  two 
solutions  are  slightly  different  analytically  but  are  in 


Figure  2.6.2  Comparison  of  the  infinite-cascade 
solution  with  existing  solutions. 
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numerical  agreement  as  seen  in  figures  2.6.3  and  2.6.4; 
the  analytical  discrepancy  could  not  be  resolved.  All  of 
the  solutions  are  in  accord  for  the  larger  phase  angles. 
Continuing  with  the  conclusions  begun  in  section  2.5.6,  it 
can  be  added  that: 

4.  For  interblade  phase  angles  larger  than  15°  the 
present  method  of  analysis  yields  closed-form  pressure- 
difference  expressions  which  are  in  agreement  with  existing 
solutions  for  the  slowly-oscillating  cascade. 

2.7  Limiting  Cases  of  the  Finite  and  Infinite-Cascade 
Solutions . 

2.7.1  The  Steady  Cascade.  All  terms  of  the  plunge  solu- 
tions are  frequency  terms;  therefore  if  k is  set  to  zero, 
the  pressure-difference  distribution  goes  to  zero  for  the 
entire  blade,  which  is,  of  course,  what  should  happen  for  a 
stationary  flat  plate  aligned  with  the  flow. 

If  k is  set  to  zero  in  the  pitching  case,  then  one 
obtains  for  the  isolated  airfoil 


(2.7.1) 


and  so  both  the  finite  (2.5.27)  and  infinite  (2.6.17) 
solutions  becomes 


(2.7.2) 


Figure  2.6.3  Comparison  of  the  infinite-cascade 
solution  with  existing  solutions. 


Figure  2.6.4  Comparison  of  the  infinite-cascade 
solution  with  existing  solutions. 


where  the  phase-angle  term  has  been  ignored  because  it  is 
meaningless  in  the  steady  case.  Equation  (2.7.2)  is  the 
classical  Ackeret  solution  [27]  for  a flat  plate  at  an  angle 
of  attack  equal  to  one  radian,  which  is  the  maximum  up-pitch 
position. 

2.7.2  The  Single-Blade  "Cascade".  If  either  cascade  is 
"squeezed"  down  to  a single  blade  by  setting  the  spacing 
parameters  m and  d to  zero,  then  the  solutions  (2.5.26,27) 
and  (2.6.16,17)  reduce  to  those  of  the  isolated  blade 

-2.Cp^  (2.7.3) 

where  the  phase  term  has  again  been  ignored. 

2.7.3  The  Infinitely-Spaced  Cascade.  The  isolated-blade 
solution  is  again  recovered  from  equations  (2.5.27)  and 
(2.6.17)  as  the  cascade  spacing  is  allowed  to  get  large. 

If  the  spacing  parameters  m and  d approach  infinity  in 
such  a way  that  the  terms  yStiM*"  and  remain  equal, 

then  equation  (2.7.3)  again  results.  The  significance  of 
these  two  terms  is  discussed  further  in  the  next  section. 

2.8  The  Unique  Inflow  Condition.  A consideration  of  the 
infinitely-spaced  cascade  in  section  2.7.3  leads  one^  to 


^The  author  is  grateful  to  Professor  M.  F.  Platzer  for 
first  asking  this  question. 

L 


ask:  For  which  cascades  does  the  relationship 


(2.7.4) 


hold?  As  seen  from  equation  (2.5.27),  for  such  a cascade 
the  pressure-difference  over  the  inflow  portion  of  the  blade 
is  merely  that  of  the  isolated  pitching  blade  times  the 
factor  C\- 

Let  <»<  be  the  Mach  angle  defined  by 


(2.7.5) 


then  it  is  not  difficult  to  show  that  for  any  cascade  having 
the  stagger  angle 


A / SlNZOt  \ 


(2.7.6) 


equation  (2.7.4)  will  also  hold.  Said  another  way,  all 
members  of  an  infinite  family  of  pitching  cascades  with  the 
stagger  angle  given  by  (2.7.6)  — and  with  no  restrictions 
on  blade  spacing  - will  have  exactly  the  same  inflow  pressure 
distribution,  i.e. 


PI: 


(2.7.7) 


This  is  termed,  by  the  author,  the  unique  inflow  condition 
for  the  slowly-oscillating  cascade.  Figure  2.8.1  depicts 


RC.U*.Cre“^) 


Figure  2.8.1  The  unique  inflow  condition. 


three  different  cascades,  each  with  the  same  stagger  angle, 
satisfying  (2.7.6).  The  pressure  distribution  for  each 
was  obtained  via  Platzer's  method-of-characteristics  pro- 
gram. The  unique  inflow  condition  is  evidenced  by  the 
exact  agreement  of  the  pressures  for  all  three  cascades 
over  the  first  portion  of  the  blade. 

2.9  Some  Comments  on  Resonance . For  a given  reference 
blade  in  an  infinite  cascade,  resonance  occurs  when  aero- 
dynamic disturbance  waves  generated  by  upstream  blade  are 
in  phase  with  the  motion  of  the  reference  blade  [12] . 
Samoylovich  [28]  shows  that  at  the  resonance  condition  the 
linear  cascade  experiences  zero  aerodynamic  damping.  The 
development  .that  follows  is  based  on  Verdon's  [12]  excellent 
explanation  of  the  resonance  phenomenon  and  provides  a means 
of  delineating  what  he  calls  the  subresonant  and  superresonant 
regions  of  the  yx-Vi  plane. 

Referring  to  figure  2.9.1,  also  taken  from  [12],  it  is 
seen  that  the  leading  edge  of  each  blade  creates  disturbances 
which  propagate  in  all  directions  at  the  freestream  speed 
of  so\ind  a^Q  . These  disturbances  are  swept  downstreaim  by 
the  freestream  flow  U , and  there  are  exactly  two  resultant 
waves  which  propagate  along  the  cascade  leading  edge  in 
the  direction  opposite  to  the  direction  of  rotation  and  at 
velocities  V,  and  V-  . From  the  law  of  cosines  one  has 


Figure  2.9.1  The  backward-traveling  waves 

causing  resonance. 
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Figure  2.9.2  The  subresonant  and  superresonant  regions. 
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The  times  required  for  a wave  from  an  upstream 

blade  (n  = -1,  -2,  ...)  to  reach  any  given  reference  blade 
(n  = 0)  are 

(2.9.2) 

The  reference  blade  lags  an  upstream  blade  by 

•to~  nss— (2.9.3) 

and  resonance  occurs  when  the  upstream  waves  are  in  phase 
with  the  reference-blade  motion,  i.e.  when 

"t-o— (2.9.4) 

or  when 

(2.9.5) 

For  a fixed  reduced  frequency  k the  superresonant  inter- 
blade phase  angles  lie  between  yn,  and  yWj , as  shown  in 
figure  2.9.2,  while  the  subresonant  phase  angles  lie  to 
either  side.  The  superresonant  region  is  important  because 
Verdon  is  of  the  opinion  that,  in  that  region,  neither  a 
finite-cascade  solution  nor  a passage-approach  infinite- 
cascade  solution  will  adequately  approximate  the  true 
infinite-cascade  solution.  In  an  effort  to  shed  some  light 

^ 


on  the  subject,  some  comparisons  are  presented  in  figures 
2.9.3  through  2.9.6.  For  cascade  A,  with  k = 0.05  , the 
superresonant  range  of  phase  angles  is  approximately  -7® 
to  -2®  , and  so  -5®  was  examined.  Figures  2.9.3  and  2.9.4 
show  that  Platzer's  finite  solution  and  the  author's  passage- 
approach  infinite  solution  are  in  agreement  with  Verdon's 
solution.  For  cascade  B,  however,  Verdon's  solution  is 
significantly  different  from  the  other  two  over  the  last 
portion  of  the  blade  for  the  superresonant  phase  angle  -3®  . 
Figures  2.9.5  and  2.9.6  show  some  high-frequency  comparisons 
for  the  phase  angle  -90®  , which  is  in  the  superresonant 
regions  of  both  cascades  A and  B.  In  these  two  cases  the 
two  solutions  — one  finite  and  one  infinite  — are  in  accord. 
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3.  THE  NONLINEAR  CASCADE 


3.1  General.  In  two  papers  Teipel  provides  the  numerical 
tools  for  the  analysis  of  low  supersonic  flow  past  a single 
oscillating  airfoil  of  finite  thickness.  In  one  paper  [23] 
he  presents  an  unsteady,  nonlinear  characteristics  method 
to  determine  the  continuous  variation  of  flow  variables. 

In  the  other  [24]  he  treats  the  oscillating  shock  waves 
present  in  the  nonlinear  problem  to  permit  determination 
of  the  discontinuous  jumps  in  the  flow  variables  across  the 
shocks  emanating  from  the  single  blade.  Using  these  tech- 
niques as  well  as  some  computational  and  theoretical  insight 
provided  by  Chadwick  [5.]‘,  who  very  ably  translated  and  applied 
Teipel 's  work,  the  present  study  is  aimed  at  finding  the 
pressures,  steady  and  unsteady,  on  the  second  blade  of  a 
parabolic-arc  cascade.  The  cascade  oscillates  in  pitch 
alone,  and  the  pressures  are  sought  over  the  entire  upper 
blade  surface  and  over  the  lower  surface  up  to  the  point  of 
impingement  of  the  trailing-edge  shock  from  the  preceding 
blade  (figure  3.2.2).  These  surface  pressures  are  computed 
by  means  of  a FORTRAN  IV  computer  program  prepared  by  the 
author  for  the  IBM- 360  computer.  The  program  listing  appears 
at  the  end  of  this  work  immediately  after  Appendix  A.  What 
follows  is  an  orderly  presentation  of  the  theory  underlying 
the  computer  program  and  the  computational  techniques  used 
in  implementing  that  theory.  Wherever  possible  the  text 
notation  will  coincide  with  that  used  in  the  program. 


3.2  Cascade  Geometry. 


a- 


P 
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3.2.1  Blade  Shape.  The  blade  shape  was  chosen  to  approx- 
imate the  blades  in  the  test  cascade  employed  by  Fleeter 
and  Riffel  [25]  who  attempted  to  measure  steady  and  unsteady 
pressures  in  a pitching  cascade.  The  airfoil  used  in  the 
computer  model  has  a flat  upper  surface  and  a lower  surface 
conforming  to  the  parabolic  arc  given  by 

A-'t  ^ lOt-V)  (3.2.1) 

where  chordwise  distance  x and  midchord  thickness  'C  are 
normalized  by  chordlength.  Figure  3.2.1  shows  the  computer 
model  for  'C  = 0.03  superimposed  on  Fleeter 's  test  blade 
and  .lists’  the  surface  coordinates  of  the  test  blade  in  per- 
cent chord.  The  greatest  discrepancy  between  the  two  shapes 
occurs  on  the  aft  portion  of  the  blade  in  the  region  behind 
the  lower-surface  shock  reflection;  however,  the  lower- 
surface  pressures  are  not  computed  in  this  region. 

3.2.2  Blade  Spacing.  The  blade  spacing  parameters  are 
input  variables  for  the  program  and  are  discussed  in  more 
detail  in  section  3.5.1.  The  physical  spacing  of  the  blades 
of  the  computer-model  cascade  — dubbed  Cascade  T — is  shown 
in  figure  3.2.2. 

3.3  Teipel*s  Method  of  Characteristics.  The  underlying 
insight  of  Teipel's  treatment  of  the  unsteady  problem  is 
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Figure  3.2.1  Test  blade  shape  and  coordinates. 


his  representation  of  the  velocity  potential  as  the  sum 


of  a steady  ‘component  and  a small-amplitude  unsteady  one 


(3.3.1) 


as  depicted  in  figure  3.3.1.  When  (3.3.1)-  is  substituted 
into  the  transonic  unsteady  small-disturbance  equation^ 


(3.3.21 


and  into  the  unsteady  boundary  condition 


^For  a discussion  of  the  limitations  imposed  by  this 
equation  see  Teipel  [23]  and  Chadwick  [5] . 


Figure  3.3.1  The  total  velocity  potential. 


where  y(x,t)  describes  the  oscillating  airfoil  surface, 
one  obtains  two  separate  problems  coupled  by  the  steady 


velocity  potential  0^  . 
and  zero-angle-of-attack 


One  problem  is  nonlinear,  steady 


(3.3.4) 


' dx 


(3.3.5) 


and  the  other  linear  and  unsteady 


O (3.3.6) 


||l  = -[.^LW<x-.)]  (3.3.7) 

where  0j^(x,y)  is  an  aunplitude  function  and  the  factor 
is  omitted.  It  is  important  to  note  that  since  the  char- 
acteristic directions  are  determined  by  the  coefficients  of 
the  highest  derivatives,  both  problems  will  have  the  same 
characteristic  directions.  The  state  variables^  for  the  two 


The  symbol  /x  is  used  both  for  the  steady  state  variable 
and  the  interblade  phase  angle,  and  it  is  normally  clear 
from  usage  which  meaning  applies. 


problems  are  essentially  the  steady  and  unsteady  velocity 
components  u^,  v^,  and 

STEADY:  M* U+O  U*  (3.3.8) 

V2  M*  (3.3.9) 

UNSTEADY:  U,,  Vj  (3.3.10) 

where 


(3.3.11) 


The  characteristic  directions  for  both  problems  are 


^ X. 


(3.3.12) 


where  the  upper  sign  refers  to  the  leffc-running  charac- 
teristics and  the  lower  to  the  right-running  character- 
istics.^ The  compatibility  relations,  however,  are  not 
the  same  for  the  two  problems . They  are 


1 2 -i- 

The  symbol  yS  is  used  both  for  the  quantity  (M  -1)^ 

and  the  right-running  characteristic  direction.  It  is 

normally  clear  from  usage  which  meaning  applies. 


where 


^0/  :s 


iCiyB  =S  JX 


(3.3.15) 


Equation  (3.3.14)  requires  some  additional  comment. 

Note  first  that  in  addition  to  the  state  variables  u^ 
and  v-^  , 0j^  is  also  present  in  the  unsteady  compatibility 
relations.  Since 


(3.3.16) 


then  must  be  found  by  integration  as  one  moves  from 

one  gridpoint  to  another.  For  the  points  A,  B,  and  C 
shown  in  figure  3.3.2,  for  example 

(3.3.17) 

S6»  + (3.3.18) 


In  performing  the  actual  computations  an  average  of  the  two 
above  expressions  is  employed. 


The  factor  in  (3.3.14)  is  straightforward,  repre- 

senting the  increment  in  7^  along  the  yS  characteristic. 
At  point  C in  figure  3.3.2  for  example  one  could  use  the 
approximation 


(3.3.19) 

However  is  a bit  more  complex,  for  while  the  increment 

dx  is  along  the  characteristic,  the  quantity 

represents  the  change  in  A with  respect  to  x along  the 
crossing  y8  characteristic  at  the  point  of  interest. 

Again  at  point  C one  could  employ 

(3.3.20) 

3 . 4 Teipel's  Treatment  of  the  Oscillating  Shock  Wave  with 
Extension  to  the  Present  Cascade. ^ 

3.4.1  The  Isentropic  Rankine-Hugoniot  Pressure  Relation . 

Consider  a normal  shoclc  moving  through  an  undisturbed 

perfect  gas  at  velocity  V relative  to  the  gas  as  shown 

s 

in  figure  3.4.1,  where  the  hat  notation  indicates  conditions 
behind  the  shock.  The  governing  continuity,  momentum  and 
energy  equations  are 


^Except  for  u,,  v.,  “u,  and  , the  velocities  appearing 
in  sections  3.4.1  through  3.4.6  are'^dimensional , in  a brief 
departure  from  the  normalization  of  all  velocities  by  the 
freestream  velocity  U . 


(3.4.1) 


yo  C- Vs’)  ^ ( Vh  - Vs  - V»,) 

P-. yiVs*  = P lVv,-Vs-V„)" 


(3.4.2) 

(3.4.3) 


Assuming  constant  specific  heats,  the  speeds  of  sound  are 
introduced,  and  the  enthalpies  are  eliminated  by 


-if  ^ 

o.  --Xj- 

(3.4.4) 

(3.4.5) 

»-! 

and,  after  a good  deal  of  manipulation,  the  change  in 
the  normal  velocities  across  the  shock  is  found  in  terms 
of  the  shock  velocity  and  sound  speed  ahead  of  the  shock 


Vv>-V»  _ 2.  

' <x  Vs 


(3.4.6) 


Equations  (3.4.1)  and  (3.4.6)  and  the  isentropic  relation 


(3.4.7) 


yield  the  Rankine-Hugoniot  pressure  relation  sought 


(3.4.8) 


r 

Figure  3.4.2  is  based  on  figure  8-3  of  [29],  and  it  shows 
that  the  isentropic  approximation  introduces  no  significant 
error  for  Fleeter 's  test  cascade  or  its  computer  model. 

3.4.2  The  Blade-Two  Shock  Position.  The  leading-edge 
shock  from  the  pitching  blade  oscillates  out  of  phase  with 
the  blade.  Let  be  the  shock  angle  of  the  steady, 

leading-edge  s’  ock  as  shown  in  figure  3.4.3.  Since  the 
shock  is  curved,  varies  with  distance  y from  the 

blade,  and  by  convention 

0<yo<^  (3.4.9) 

The  shock  angle  at  any  point  on  the  oscillating  wave  can 
be  written  as  a small  variation  from  the  steady  one^ 

1 

y a 'jfo (3.4.10) 

i 

I In  order  to  locate  the  oscillating  shock  relative  to  the 

steady  one,  Teipel  introduces  the  complex  displacement 
function  g(y)  , also  shown  in  figure  3.4.3,  and  so  the  x 
coordinate  of  a shock  point  is  approximately 

COT'tfo  + (3.4.11) 


The  symbol  gamma  used  here  and  in  figures  3.4.4  and 
3.4.5  represents  the  total  shock  angle.  All  other  uses  of 
the  unsubscripted , unsuperscripted  gamma  stand  for  the 
dimensionless  ratio  of  specific  heats. 


r 


where  steady  coordinates  of  the  blade's 

leading  edge.  Taking  the  total  differential  of  (3.4.11) 
and  noting  that  along  the  shock 


ax 


COT  'H 


(3.4.12) 


results  in 


s (3.4.13) 

thus  relating  shock  angle  and  shock  displacement. 

3.4.3  The  Velocity  Components  Upstream  of  the  Blade-Two 
Shock  Wave . Ahead  of  the  shock  the  velocity  components 
have  both  steady  and  unsteady  parts 

w.*  = "U  -V  Uo  -f  u,  (3.4.14) 

V*  = V,  T V,  (3.4.15) 

where  the  barred  velocities  are  dimensional  counterparts  of 
the  dimensionless  quantities  u^,  v^,  u^^  and  v^  already 
introduced  in  equations  (3.3.8,9,10).  However,  the  flat  upper 
surfaces  of  the  blades  of  cascade  T are  aligned  with  the  flow 
in  the  steady  position,  and  so,  upstream  of  the  shock 

Uo  = "Vo  ST  O 


I 


(3.4.16) 


Referring  to  figure  3.4.4  one  has  for  the  normal  and 
tangential  velocity  components  ahead  of  the  shock 

= + U cos  V*  v'  -»•  (u,  + V,  COS  jfo'i  (3.4.17) 

V^  = lTcos>f,-XJ  s\N>5,<'  -►  (u,  cos^o-  V,  S»N  (3.4.18) 

3.4.4  The  Blade-Two  Shock-Motion- Induced  Velocity. 
Relative  to  the  gas  ahead  of  it  the  shock  is  moving  at 
the  velocity 


V4=  -V„-»rV' 


(3.4.19) 


where  V is  a normal  velocity  induced  by  the  motion  of 
the  shock,  and  is  seen  from  figure  3.4.5  to  be 


V'=|^  2>\N-» 

O TL 


(3.4.20) 


or,  from  (3.4.11) 


V'  = UK  ^ e.''-''*-  (3.4.21) 

3.4.5  The  Velocity  Jump  Conditions  Across  the  Blade-Two 
Shock.  Using  equation  (3.4.13)  to  eliminate  y'  from  the 
expressions  for  the  upstream  velocity  components  (3.4.17,18), 
they  become 


8 


Figure  3.4.5  The  shock-motion-induced  velocity. 


V* 


sUsiKI  Vo  ■*■(■□  COi^o  U,  SIN){o  + V,COS1Jo)e'’•''*■ 

a^ 


:'U  COS«o+C-'US»N^Vo'^  •*-  u,  cos  Vo  - V,  SlNW*) 


(3.4.22) 

(3.4.23) 


Note  that  each  of  the  above  velocities  has  three  component 
parts.  The  first  term  in  each  expression  arises  from  the 
steady  shock  position,  the  second  arises  from  the  oscillatory 
motion  of  the  shock  about  that  steady  position,  and  the 
third  and  fourth  terms  arise  from  the  unsteady  disturbances 
ahead  of  the  shock  which  are  generated  by  the  first  blade's 


motion. 


The  components  behind  the  shock  are  found  by  noting 
that  the  tangential  component  does  not  change 


(3.4.24) 


and  by  using  (3.4.6,19,21,22)  to  obtain 


“SX.  S1N*X.  ll- ^ (3.4.25) 

+ V,  cos 


where 


Mh=M  sin  Vo 


(3.4.26) 


It  is  convenient  to  once  again  return  to  the  dimensionless 
x-y  velocity  components  u'^^  and  , obtaining  from  the 
unsteady  terms  of  (3.4.24,25) 


. . 

u,  = *^'5^  V*  *^2.  ^ W.I  -»•  vn.^  V, 


v»  = hj  u,  -► 


(3.4.27) 


(3.4,28) 
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Ift»  = — Ceos  i-a-o  S.IN* 

»+ 1 Mn 


W,s  -2:—  »^CV-V-  — V'  COT-lCo  SlN*-y» 
y-v*  Mn 


lrTj=  Yn^. 


(3.4.33) 


(3.4.34) 


(3.4.35) 


(3.4.36) 


1 


Note  that  the  above  coefficients  are  functions  of  the 
local  shock  angle  y ^ and  so  are  not  constant  along  the 
shock . 

If  no  unsteady  disturbances  exist  upstream  of  the  shock 
(Uj  = = 0)  then  (3.4.27,28)  reduce  to  Teipel's  single- 

blade results.  If,  in  addition,  the  shock  does  not  oscillate 
(g  = dg/dy  = 0)  , then  'u'^^  and  are  identically  zero 

indicating  that  there  are  no  unsteady  disturbances  behind 
the  shock.  Equations  (3.4.27,28)  represent  an  extension  of 
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Teipel's  work  from  the  single  blade  to  the  cascade,  and 
they  were  first  obtained  by  Chadwick  for  a cascade  of 
wedges.  It  is  noted  for  completeness  that  the  above 
equations  would  apply  equally  well  to  an  upper-surface 
shock  wave  if  the  signs  of  m^^,  m^ , and  n^  were 

reversed. 


3.4.6  The  Pressure  Jump  Conditions  Across  the  Blade-Two 
Shock.  An  additional  extension  of  Teipel's  work  was  accom- 
plished by  the  author  in  order  to  develop  a more  accurate 
method  of  determining  the  flow-variable  discontinuities 
than  that  employed  by  either  Chadwick  or  Teipel.  The 
improved  method  is  treated  in  more  detail  in  section  3.8, 
but  it  is  based  on  the  extension  of  the  Rankine-Hugoniot 
pressure  relation  (3.4.8)  to  the  cascade  via  equations 
(3.4.19,20,21),  resulting  in 


(3.4.37) 


where 


jb.s  mJ:  S>N 


(3.4.38) 


4-V  ..a. 


= -rr-r  Mm 


(3.4.39) 


(3.4.40) 
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MS  COT  )S„ 


(3.4.41) 


Equation  (3.4.37)  relates  the  pressure  jump  across  the 
shock  to  the  shock  displacement  and  the  upstream  unsteady 


velocities.  In  the  single-blade  case  u^^  and  are 


zero,  and  Teipel's  result  is  recovered. 


3.4.7  The  Boundary  Conditions  for  the  Shock  Equations. 
In  section  3.8,  the  flow  variables  behind  the  shock  are 
obtained  by  the  simultaneous  solution  of  the  difference 
equations  resulting  from  (3.3.14)  and  the  shock  equations 
(3.4.27,28,37).  In  order  to  start  the  solution  process, 
then,  several  quantities  must  initially  be  known  at  the 
poin.t  of  intersection  of  the  shock  and  blade,,  i.e.  at 


(x^,y^) . From  figure  3.4.3  it  is  clear  that 


(Wo-)  s-O 


(3.4.42) 


Also,  immediately  behind  the  shock  the  flow  tangency 
condition  for  the  second  blade  must  hold 


V, 


* = - [l  -►  A-t  (3.4.43) 


With  Vj^  known  equation  (3.4.28)  provides 


^ a CVi-V)jU,  - A-t  (7c./j;o) 

<7 


(3.4.44) 
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and  (3.4.44)  used  with  (3.4.27,37)  gives  the  remaining 
boundary  conditions 


I 

[ 


-*• = 


aJu 


(3.4.45) 


K>- 


R'V  V. 

Wl 


CXo^ 


(3.4.46) 


Of  course  for  equations  (3.4.43)  through  (3.4.46)  the  hatted 
quantities  are  just  behind  the  shock  at  ' 

Equation  (3.4.45)  was  derived  by  Chadwick  while  (3.4.46) 
is  the  author's  result. 

! 

3.4.8  The  Velocity  Potential  Jump  Conditions  Across  the 
Blade-Two  Shock.  Landahl  [26]  shows  that  the  total  velocity 
potential  (3.3.1)  is  continuous  across  the  leading-edge 
shock,  but  this  does  not  imply  that  either  the  steady  or 
unsteady  component  is  individually  continuous . Teipel  and 
Chadwick  choose  to  assume  that,  at  the  first  shock  point 

S&,C-Xo,^o)=:  (3.4.47) 

The  author,  on  the  other  hand,  chooses  to  find  the  initial 
velocity  potential  jump  at  the  point  ^ 

perfect  gas  with  constant  specific  heats,  it  can  be  shown  that 


1 


90 


01- j ) 


(3.4.48) 


/V 

and  since  and  P/Poo  known  at  from 


equations  (3.4.45,46),  then  so  is  0^  . Figure  3.4.6  is  a 
plot  of  computer  data  for  cascade  T,  and  it  shows  the 
discontinuity  in  0^^  along  the  shock  as  a function  of 
distance  from  the  blade.  Note  that  the  initial  discon- 
tinuity is  small  justifying  Teipel's  assumption  (3.4.47), 
but  equation  (3.4.48)  eliminates  having  to  make  that 
assvimption  at  all. 


3.5  Constructing  the  Characteristic  Grid.  Because  the 
characteristic  directions  are  the  same  for  both  the  steady 
and  unsteady  problem,  the  characteristic  grid  can  be  con- 
structed and  employed  to  find  both  the  steady  and  unsteady 
velocities.  The  grids  above  blade  one  (zone  A)  and  above 
blade  two  (zone  C)  are  composed  of  two  families  of  straight, 
parallel  characteristics.  The  mesh  below  blade  two  (zone  B) 
however,  is  composed  of  one  straight,  divergent  family  (yS) 
and  one  curved  family  («(). 


3.5.1  The  Straight  Characteristic  Grids  Above  Blades  One 
and  Two.  Since  the  upper  blade  surfaces  generate  no  steady 
disturbances,  it  is  seen  from  (3.3.8)  that,  in  zones  A and  C 


(3.5.1) 


9] 


Figure  3.4.6  The  discontinuous  jump  in  the  unsteady 
velocity  potential  across  the  leading- 
edge  shock. 
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and  so  the  characteristics  in  those  zones  are  simply  Mach 
waves  (figure  3.5.1).  Each  blade  has  21  equispaced  upper- 
surface  gridpoints , and  the  program  requires  that  the  leading 
edge  of  the  second  blade  be  positioned  at  one  of  the  straight- 
mesh  gridpoints.  This  is  accomplished  via  the  input  parameters 
K and  N which  fix  the  cascade  blade  spacing  according  to 


Y 


B = O.OiS 

r 


(3.5.2) 


L 


XB  = VS  yB  O.OS- M (3.5.3) 

For  cascade  T,  K and  N are  16  and  7 respectively. 
Figure  3.5.1  also  shows  the  characteristic  numbering  conven- 
tion employed  in  the  program,  so  that  a gridppint  may  be 
located  by  two  characteristic  coordinates 

(3.5.4) 

Note  that  the  straight  mesh  is  computed  as  if  both  blades 
were  flat  plates.  The  x-y  coordinates  of  the  straight 
mesh  are  referred  to  as  XI,  Yl  in  the  computer  program. 

3.5.2  The  Curved  Characteristic  Grid  Below  Blade  Two. 

In  zone  B (figure  3.5.2)  the  characteristic  grid  is  non- 
linear and  must  be  constructed  iteratively.  Figure  3.5.2 
reveals  the  numbering  convention  employed  for  the  curved 
mesh.  The  curved-mesh  x-y  coordinates  are  referred  to  as 
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Figure  3.5.1  The  straight  characteristic  grid. 


X2,y2  in  the  program  to  distinguish  them  from  those  of 
the  straight  mesh.  Before  proceeding  with  the  grid  con- 
sturction,  some  useful  relations  are  obtained. 

For  an  o<  characteristic  passing  through  the  shock 
from  zone  A to  zone  B (figure  3.5.3) 


/'a  ' 


(3.5.5) 


and  since 


A 


IS  zero 


(3.5.6) 


Now  consider  two  points  on  the  same  ^ characteristic, 

also  illustrated  in  figure  3.5.3.  From  the  c<  compatibility 

relation 


?v(X^3-)  -yw  (X,T)  = 
'ACX^K)^-yuCX,K)  ^ 


and  from  the  yS  compatibility  relation 


and  so  it  must  be  that 


>CX,T)  = A K)  JA  (X,X)  = yx  LX,K  ) 


(3.5.7) 


(3.5.8) 


A (X,  y* (X,<J)  = A (X,K)^»  ^ tX,K)  (3.5.9) 


(3.5.10) 
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which  proves,  in  addition,  that  the  yS  family  is  composed 
of  straight  lines. 

The  steady  shock  angle  at  any  point  along  the  shock  is 
found  by  averaging  the  angles  of  the  characteristic  lines 
before  and  after  the  shock,  i.e. 

* T t (3.5.11) 

Taking  the  tangent  of  (3.5.11)  and  applying  the  sum  and 
double-angle  trigonometric  formulae  leads  to  a quadratic 
equation  in  tan 

TAN*-J(^  -v  Z(V^*'^’'‘-0.TAN  =lO  (3.5.12) 

The  solution  of  (3.5.12)  is  rather  cumbersome  and  is 
adquately  approximated  by  the  simpler  expression 

TAN  ■JJ*-  Vi  (3.5.13) 

In  the  computational  implementation  of  (3.5.13),  referring 
to  figure  3.5.3  again,  the  characteristic  eingles  at  adjacent 
points  behind  the  shock  are  first  averaged  together,  and 
the  result  is  averaged  with  the  zone-A  angle,  i.e. 


TAN  -Jj*  = 


(3.5.14) 


3. 5. 2.1  Computing  the  Point  (1,1).  The  X2-Y2  coord- 


inates of  (1,1)  are  of  course  those  of  the  blade-two 
leading  edge.  yU(l,l)  is  found  from  the  steady  boundary 
condition  (3.3.5)  and  equation  (3.3.9)  while  7v.(l,l)  follows 
from  (3.5.6). 

3 . 5 . 2 . 2 Computing  the  Points  (1,2)  and  (2,2).  The 
points  (1,2)  and  (2,2)  are  established  iteratively  as 
follows.  Program  subroutines  employed  are  indicated  in 
parentheses,  and  figure  3.5.4  is  provided  for  reference. 

1.  Establish  the  reference  point  XA  on  the  blade, 
close  to  the  leading  edge 

XA=  O.OOS 

2.  Compute  yU  at  XA  , which  is  equal  to  yM(l,2), 
using  the  steady  boundary  condition  (3.3.5)  and  (3.3.9). 

3.  Compute  X at  XA  , which  is  equal  to  X(l,2), 
using  equation  (3.5.6). 

4.  (Subroutine  SOLVE)  Locate  (1,2)  at  the  intersection 
of  a line  through  (XA,  Y2(l,l))  with  slope 

-> 

and  a line  through  (X2(l,l),  Y2(l,l))  with  slope 

2± 


k 
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Figure  3.5.4  Computing  the  shock-front  points. 


I 5.  (Subroutine  ITER)  Guess  at  the  position  of  (2,2) 

i 

; on  the  blade 

j = X2.tv0 

4 

I 

i 

I, 

6.  (Subroutine  ITER)  Increment  that  guess 
71(2,2)  = X2(2,2')  -1- 0.00\ 

7.  (Subroutine  ITER)  Compute  yU(2,2)  from  (3. 3. 5, 9). 

8.  (Subroutine  ITER)  Compute  ^(2,2)  from  (3.5.6). 

9.  (Subroutine  ITER)  Solve  for  X2(2,2)  from  the 
slope  of  the  ^ characteristic  through  (1,2)  and  (2,2) 

X2(2,2)~X2t>,Z)  _ Z 

10.  (Subroutine  ITER)  Does  X2(2,2)  = X2(2,2)? 

YES:  GO  TO  11;  NO:  GO  TO  6. 

11.  STOP. 

3 . 5 . 2 . 3 Computing  Subsequent  Shock  Front  Points. 
Subsequent  points  (1,1)  along  the  shock  are  also  found 
iteratively.  The  technique  follows,  and  figure  3.5.4 
applies. 

1.  Relocate  XA  on  the  blade 

XAv-s  X/k  -V-  DXA. 


inn 


r 


where  DXA  is  based  on  the  first  position  of  XA 


DXA=: 


XZ  (Z.Z)  - XA 

^ 15 


2.  (Subroutine  SIIOKPT)  Guess  ^(1,1) 

= A(l;X-0 

3.  (Subroutine  SHOKPT)  Increment  that  guess. 

ACt,X)  = X(l,X^  -»•  o.oo» 

4.  (Subroutine  SHOKPT)  Compute  using  and 

X at  (2,1)  , which  are  the  same  respectively  as  yU(2,2) 
and  7V(2,2) 

5.  (Subroutine  SOLVE)  Locate  (1,1)  at  the  intersection 
of  a line  through  (XA,  Y2(l,l))  with  slope 


_zL 


x(i,xV'^ 


and  a line  through  (X2(1,I-1),  Y2(1,I-1))  with  slope 


— 4- 
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6.  (Subroutine  SHOKPT)  Compute  yu(l,I)  from  the 
steady  boundary  condition  (3.3.5)  at  XA  and  equation 
(3.3.9) . 


7.  (Subroutine  SHOKPT) 
YES:  GO  TO  8;  NO;  GO  TO  3. 

8.  STOP. 


= lAil,!) 
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3 . 5 . 2 . 4 Computing  the  Curved-Mesh  Interior  Points.  A 
point  (I,J)  interior  to  the  curved  mesh  is  located  (figure 
3.5.5)  at  the  intersection  of  a line  passing  through 
(I-1,J)  with  slope 


\ 

and  a line  through  '(1,1)  with  slojie 


3 . 5 . 2 . 5 Computing  Subsequent  Blade-Surface  Points. 
Subsequent  points  of  the  form  (1,1)  on  the  lower  surface 
are  computed  in  a manner  exactly  analagous  to  the  computation 
of  the  point  (2,2)  as  described  in  steps  5 through  11  of 
section  3. 5. 2. 2. 

3.6  Computing  the  Steady  Velocities  and  Pressures.  With 
yU  and  % known  throughout  the  grid,  the  steady  velocities 
follow  directly  from  (3. 3. 8, 9).  The  steady  pressures  on 


Figure  3.5.5  Computing  the  curved-mesh 
interior  points. 


the  upper  surface  of  blade  two  are  zero  and  on  the  lower 
surface  are  given  at  the  surface  points  (1,1)  by 

s - 2.  Uo (3.6.1) 

3.7  Computing  the  Unsteady  Flow  Variables  at  the 
Straiaht-Mesh  Grid  Points. 

3.7.1  The  Blade-One  Upper-Surface  Grid  Points.  In  the 
straight  mesh  the  unsteady  velocities  and  velocity  potential 
are  designated  by 

uj^  = UllR  ♦ i mil  vj^  a VllR  + i VllI 

= FlR  -*■  i PlI  (3.7.1) 

- J 


r 

\ 

The  unsteady  boundary  condition  for  blade  one  (3.3.7) 
allows  determination  of  the  velocities 

VllR(I,I)  = -1  V11I(I,I)  = -k[xi(I,I)  - b]  (3.7.2) 

Using  the  yg  compatibility  relation  (3.3.14)  in  difference 
form  results  in  a system  of  linear  equations  for  the 
unknowns  U11R(I,I),  U11I(I,I)  . Figure  3.7.1  A applies, 

and  the  equations  appear  in  matrix  form  in  section  A.l  of 
the  appendix.  The  velocity  potential  is  found  by  integration 
of  the  velocities,  using  equation  (3.3.18);  see  section  A. 2 
of  appendix  A. 

3.7.2  The  Straight-Mesh  Field  Points.  As  seen  in  figure 

. 3.7.1  B,  grid  points  not  on  the  blade  require  use  of  both 

compatibility  relations.  The  equations  resulting  for  the 
unknown  u^^  and  v^^  velocities  at  each  field  point  (I,J) 
are  given  in  A. 3 of  the  appendix.  Equations  (3.3.17,18) 
are  averaged  to  find  the  velocity  potential;  see  section 
A. 4 of  the  appendix. 

3.7.3  The  Grid  Points  Immediately  Upstream  of  the  Shock. 
Figure  3.7.1  C shows  that  the  shock  points  lie  in  the  straight 
mesh  and  therefore  the  flow  variables  at  each  shock  point  — 
upstream  of  the  shock  — are  found  by  means  of  an  interpolation 
algorithm.  The  distances  RI  in  figure  3.7.1  C are  computed 
and  any  desired  quantity  Q(1,K)  is  found  from  the  values 

of  Q at  the  four  surrounding  points 


^ f\  A 


Figure  3.7.1  Computing  the  straight-mesh  and 
curved-mesh  grid  points. 
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(3.7.5) 


_ RR>  QCI-t.7-0-*-RRZ  Q(X.-;r^->-KR4-  QCX.7-rt 

' ” RD 


where 


RR^s:yfR.X  RD  = 5II  (3.7.4) 

1^7 

3. 8 Determining  the  Jumps  in  the  Flow  Variables  Across  the 
Blade-Two  Shock.  It  is  in  the  determination  of  the  jump 
conditions  that  the  author  employs  a technique  — alluded  to 
in  section  3.4.8  — that  is  different  from  the  method  used 
by  Teipel  in  the  single-blade  case.  The  essence  of  Teipel's 
approach  is  to: 

1.  Assume  that  the  velocity  potential  does  not  change 
across  the  shock  at  the  first  point  (3.4.47). 

2.  Use  the  shock  velocity  equations  (3.4.27,28),  with 
boundary  conditions  (3.4.42,43,45),  in  conjunction  with  the 
yfi  compatibility  relation  (3.3.14),  to  solve  for 

(behind  the  shock)  and  g at  each  shock  point. 

3.  Use  equation  (3.3.18)  behind  the  shock  to  recover 
0^  at  each  shock  point. 

In  contrast,  the  method  employed  in  the  present  study  is  to: 

1.  Use  (3.4.46,48)  to  determine  0^  behind  the  shock 
at  the  first  shock  point. 
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2.  Use  the  shock  velocity  equations  (3.4.27,28),  the 
shock  pressure  equation  (3.4.37),  the  yS  compatibility 
relation  (3.3.14)  and  the  boundary  conditions  (3.4.42,43, 

45,46)  to  solve  for  “u^,  g and  ^/P  at  each  shock  point. 

3.  Use  equation  (3.4.48)  to  find  0^  from  the  pressure 
jump  at  each  shock  point. 

The  latter  method  results  in  a larger  system  of  equations 
at  each  point,  but  no  assumptions  concerning  the  continuity 
of  the  velocity  potential  are  required. 

3.8.1  Details  of  the  Jump-Condition  Computations. 

Step  2 above  requires  some  additional  comment.  The  shock 
equations  (3.4.27,28,37)  are  applied  in  difference  form  from 
one  shock  point  to  the  next.  They  are  the  last  six- equations 
in  the  ^aatrix  equation  given  in  appendix  section  A. 5.  The 
y8  relation  supplies  the  first  two  equations  in  that  matrix 
equation.  Note  from  figure  3.7.1  D that  using  the  J3  relation 
at  shock  point  (1,1)  requires  two  tacit  assumptions.  First, 
the  flow  variables  at  the  point  where  the  yg  characteristic 
through  (1,1)  leaves  the  blade  are  essentially  those  at 
the  point  (1,1)  immediately  behind  the  shock.  Second,  the 
flow  variables  on  the  y8  characteristic  through  (1,1+1) 
and  just  behind  (1,1)  are  about  the  same  as  those  at 
(1,1)  behind  the  shock. 

In  the  computer  program  the  flow  variables  at  shock 
point  (1,1)  , ahead  of  the  shock,  have  the  form 


= U12(1,I)  - i Ul2I(l,i) 


(3.8.1) 


while,  at  the  same  point  immediately  behind  the  shock, 
they  are  designated 

Ml  a URJ(I)  + i UIJ(I)  (3.8.2) 


until  stored  in  the  form  of  (3.8.1). 


3 . 9 Computing  the  Unsteady  Flow  Variables  at  the  Curved- 
Mesh  Grid  Points. 

3.9.1  The  Blade-Two  Lower-Surface  Grid  Points.  The 
unsteady  boundary  condition  for  the  second  blade  provides 
the  lower-surface  velocities  (figure  3.9.1  A).  The 

equations  are  given  in  section  A. 6 of  the  appendix.  The 
Uj^  velocities  are  found  via  the  oi  compatibility  relation 
(section  A. 7)  and  the  velocity  potential  follows  from 
equation  (3.3.17)  (appendix  section  A. 8), 


3.9.2  The  Curved-Mesh  Field  Points.  Both  compatibility 
relations  are  used  at  a field  point,  as  seen  from  figure 
3.9.1  B.  The  resultant  matrix  equation  appears  in  appendix 
section  A. 9.  The  velocity  potential  is  found  from  equations 
(3.3.17,18)  (appendix  section  A. 10). 

With  the  curved  mesh  completely  computed  the  unsteady 
pressures  at  blade-two,  lower-surface  grid  points  are  given 
by 


= -2[ui2R(I,I)  - k F2I(I,I)] 
-2i[ui2I(I,I)  ♦ k F2R(I,I)] 


(3.9.1) 
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Figure  3.9.1  Computing  the  curved-mesh  grid  points. 


3.10  Results  and  Conclusions. 

3.10.1  The  Isolated-Blade  Case.  If  the  unsteady  boundary 
condition  along  the  first  blade  is  set  to  zero,  the  second 
blade  sees  freestream  conditions  and  its  unsteady  pressures 
can  be  compared  with  Teipel's  isolated-blade  results  [23]. 
This  is.  accomplished  in  figure  3.10.1,  where  it  is  noted 
that  a small  discrepancy  between  the  two  solutions  exists. 
This  discrepancy  is  attributable  in  part  to  the  different 
techniques  used  by  Teipel  and  the  author  in  solving  the 
shock  equations  (see  section  3.8).  But  it  is  also  noted 
that  Chadwick  [5]  too  encountered  a discrepancy  when  he 
applied  Teipel's  method  to  the  single,  biconvex  airfoil. 
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Teipel  points  out  that , for  the  single  blade,  as  k 
gets  large  the  linear  and  nonlinear  solutions  approach 
each  other.  This  phenomenon  is  also  evident  for  the  second 
blade  in  cascade  T,  as  seen  in  figure  3.10.2  which  compares 
Platser's  linear  solution  [17]  with  the  nonlinear  results  of 
this  investigation  for  a k of  1.5  . 


3.10.2  The  Steady  Pressure  Results.  Figure  3.10.3  shows 
the  steady  pressures  for  the  lower  surface  of  blade  two. 

It  is  seen  that  agreement  between  the  experimental  results 


Ob. 

u 


Figure  3.10.2  High-k  comparison  of  the  linear 
and  non-linear  solutions. 

of  Fleeter  [25]  and  the  numerical  results  of  the  computer 
model  is  good  at  midchord  and  fair  at  the  leading  edge. 

The  program  results  compare  favorably  with  theoretical 
predictions  from  second-order  Prandtl-Busemann  theory  and 
I a third-order  theory  due  to  Ferri  [30]  which  takes  entropy 

losses  into  account.  This  latter  comparison  again  demon- 
strates that  the  computer  solution  is  not  limited  by  its 
I use  of  the  isentropic  Rankine-Hugoniot  pressure  relation 

(3.4.8). 

3.10.3  The  Unsteady  Pressure  Results.  Figures  3.10.4 
through  3.10.12  show  the  computer  results  of  the  present 
study,  the  experimental  data  from  Fleeter 's  investigation, 
and  the  linear  results  of  Platzer's  method-of-characteristics 
I program  for  nine  interblade  phase  angles.  The  pressure 

k.  


1 


Figure  3.10.3  Blade-two  lower-surface 
steady  pressures. 


data  are  presented  as  an  amplitude  and  phase  angle  relative 
to  the  motion  of  the  second  blade.  The  phase  angle  is 
measured  from  the  point  where  the  second  blade  attains 
maximum  up  pitch  to  the  point  of  peak  blade  pressure.  Along 
the  upper  surface  of  blade  two  the  linear  and  nonlinear 
numerical  results  agree  exactly  — because  the  upper  surface 
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is  flat  — providing  an  additional  check  of  the  nonlinear 
program.  Along  the  lower  surface  the  nonlinear  effects 
show  up  as  differences  between  the  two  predicted  pressure 
distributions  — differences  which  vary  with  interblade  phase 
angle.  The  predicted  nonlinear  pressure  autiplitude  is,  in 
all  cases  examined,  at  least  as  large  as  the  linear  predic- 
tions, which  is  as  expected.  In  figure  3.10.5  (^=  -4.4“) 
the  linear  and  nonlinear  amplitude  predictions  are  nearly 
coincident,  while  in  figure  3.10.9  130.6“)  they  differ 

considerably.  Figure  3.10.13  is  a plot  of  the  difference 
between  the  nonlinear  and  linear  amplitude  predictions  versus 
the  turning  angles  experienced  by  the  flow  at  the  leading 
edge  of  blade  two.  It  shows  that  the  nonlinear  computer 
predictions  are  consistent  since  the  variation  from  the 
linear  solution  increases  with  the  flow  turning  angle. 

That  turning  angle  is  given  by 


T.A^.=  57.  3*  COS/*) 


(3.10.1) 


when  blade  two  is  in  the  maximum  up-pitch  position.  Note 
that  the  large  turning  angles  in  figure  3.10.13  result  from 
the  choice  of  a one  radian  maximum  pitch  amplitude  for 
nondimensionalization  purposes. 

Comparison  of  the  nonlinear  computer  predictions  with 
Fleeter 's  experimental  data,  however,  reveals  that  the 
theory  and  the  experimental  results  do  not  consistently  agree. 
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Percent  Chord 

Figure  3.10.4  Blade-two  unsteady  pressures. 
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Figure  3.10.7  Blade-two  unsteady  pressures. 
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Figure  3.10.11  Blade-two  unsteady  pressures. 
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Interblade  phase  angle  (degrees) 

Figure  3.10.13  Nonlinear  variations  from  linear  theory 

In  figure  3.10.11,  for  example,  there  is  adequate  upper- 
surface  agreement  with  the  theory,  but  the  lower-surface 
pressure-amplitude  comparison  is  poor.  In  fact,  in  no  case 
is  there  agreement  between  data  and  theory  for  both  phase 
and  amplitude  along  the  lower  surface.  The  lower-surface 
experimental  data  appear  to  be  somewhat  inconsistent.  They 
sometimes  vary  widely  from  one  point  to  the  next  for  the 
same  interblade  phase  angle  (figure  3.10.6)  and  they  also 
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vary  a great  deal  for  rather  small  changes  in  interblade 
phase  angle  (figures  3.10.10,11,12).  In  addition,  in  each 
of  the  figures  there  is  at  least  one  amplitude  data  point 
for  the  nonlinear  (lower)  surface  that  is  below  the  corres- 
ponding linear-surface  data  point.  This  phenomenon  is 
contrary  to  expectation  because  the  thickness  surface  must 
generate  a larger  magnitude  disturbance  than  the  flat  surface. 
Certainly,  additional  experimental  data  are  required  to 
adequately  evaluate  the  theoretical  methods. 

3.10.4  Conclusions.  In  summary,  it  is  concluded  that  the 
extension  of  Teipel's  single-blade  methods  to  the  nonlinear 
cascade  results  in; 

1.  Steady  pressure  predictions  which  are  in  accord 
with  other  nonlinear  theories  and  which  are  supported,  to 
some  extent,  by  experimental  results. 

2.  Unsteady  pressure  predictions  which  exhibit  expected 
trends  relative  to  linear  theory,  but  which  are  not  consis- 
tently in  agreement  with  experimental  data  presently  available. 
The  significant  influence  of  thickness  on  the  unsteady 
pressures  is  apparent  in  both  the  theory  and  data. 

3.  A numerical  approach  to  the  unsteady  cascade  which 
appears  to  have  the  potential  to  treat  the  effects  of  blade 
thickness,  blade  camber,  and  steady  angle-of-attack  loading 
to  provide  a more  realistic  model  of  the  fan  and  compressor 
blades  in  a turbomachine. 
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APPENDIX  A 


A,1  Equations  for  the  Unknown  ui  Velocity  at  Blade-One 


Surface  Points. 

1.  References:  section  3,7.1;  figure  3.7.1  A. 

2,  Auxiliary  constants: 

Al  = 2kM2(0.025)A.  A2  = k2M2(0.025)/X. 


A3  = !•*•  0.0125  A2 
A5  = 1-  0.0125  A2 
3,  Matrix  equation: 


A4  = A3/7s5[ 
A6  = 


Al/2 


-Al/JI  [B.1R(I,I)] 
A5^  mid,  I) 


A3  U11R(  1-1,1)  -h 
A3  U11I(I-1,I)  - 


♦ Al  UllI(I-l,I)/2  - A4  V11R(I,I)-¥A6  VUR(I-1,I)  + 
j^Al  UllR(I-l,I)/2  - A4  vni(I,I)-»-A6  VllI(I-l,I)*»- 

+ A2  F1R(I-1,I) 

-►A2  F1I(I-1,I) 


Ao2  Equations  for  f6^  at  Blade-One  Surface  Points. 


1.  Reference:  section  3.7,1;  figure  3.7.1  A. 

2,  Equations: 

F1R(I,I)  = F1R(I-1,I)  -vlullRd,!)  UllRCl-l,!)  - 

- [vllRCI,!)  ’¥  VllRd-1,1  (0.0125) 

F1I(I,I)  = FlId-1,1)  {ullI(I,I)  + UllId-1,1)  - 

- [yiii(i,i)  -fviiid-i,!  (0.0125) 
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A, 3 Equations  for  the  Velocities  ui  and  vi  at  Straight -Mesh 
Grid  Points. 


1.  References:  section  3,7,2;  figure  3.7,1  B, 

2,  Auxiliary  constants: 


A1  = 2kM2(0.025)/7v 
A3  = 1 - 0.0125  A2 
A5  = 

3,  Matrix  equation: 


A2  = k2M2( 0.025)/^ 
A4  = 1 ^ 0.0125  A2 

A6  = A4/>.Vi 

A 


"*A3 

-Al/2 

-A6 

Al/2 

A3 

0 

A3 

-Al/2 

A6 

Al/2 

A3 

0 

o"” 

‘1jiir(i,jT 

A6 

U11I(I,J) 

0 

V11R(I,J) 

A6 

V11I(I,J) 

A4  U11R(I,J-1)  + 
A4  U11I(I,J-1)  - 
A4  U11R(I-1,J)  ^ 
A4  U11I(I-1,J)  • 


•»»  Al  U11I(1,J-1)  - A5  V11R(I,J-1)  A2  FlR(l,J-l) 

2. 

- U11R(I,J-1)  - A5  V11I(I,J-1)  -V  A2  FlI(I,J-l) 
UllI(I-l,J)  -V  A5  V11R(I-1,J)  -¥  A2  F1R(I-1,J) 

2 

- ^ U11R(I-1,J)  -*•  A5  V11I(I-1,  J)  A2  F1I(I-1,  J) 

Z — 


A, 4 Equations  for  at  Straight -Mesh  Grid  Points. 

F1R(I,J)  « i[FlR(I,J-l)  + FlRCl-l.jJ  0,0125  UllR(I,J) 

-f  0.00625  |ullR( I,  J-1)  -t  U11R(I-1,J^ 

+ 0.00625  lyllRC I,  J-1)  - V11R(I-1, 

FlI(I,J)  = i[FlI(I,J-l)  F1I(I-1,j3  ^ 0.0125  UllI(I,J)  ^ 
0.00625  |irLlI(I,J-l)  + U11I(I-1,  J)]  -► 

0.00625  [vllI(I,J-l)  - V11I(I-1, 

1,  References:  section  3.7.2;  figure  3,7,1  B, 
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A.5  Equations  for  the  Jmp  Conditions  across  the  Blade-Two 
Leading-Edge  Shock. 

1.  References:  section  3.8.1;  figure  3.7.1  D. 

2.  Atixiliary  constants: 

A1  = 2/  [^(1,1)  +i^Cl,I-liI 

A2  = 4kM2[p(l,I)  - X2(1,I-1)]/  [7^(1, 1)  +■^(l,I-l)] 
A3  = A2  k/2 

DY  = [y2(1,I)  - Y2(1,I-i3/2 

3.  Matrix  equations: 


1 

-A2/2 

Al 

0 

-A3/2 

0 

0 

0 

URJ(I) 

A2/2 

1 

0 

Al 

0 

-A3/2 

0 

0 

UIJ(I) 

1 

0 

0 

0 

0 

0 

-mi/DY 

m2 

VRJ(I) 

0 

1 

0 

0 

0 

o' 

-m^/DY 

V.IJ(I) 

0 

0 

1 

0 

0 

0 

-ni/DY 

"2 

FRJ(I) 

0 

0 

0 

1 

0 

0 

“"2 

-nj^/DY 

FIJ(I) 

-tM^ 

0 

0 

0 

0 

VM^k 

-Pl/DY 

P2 

GR(I) 

0 

-»m2 

0 

0 

-^M^k 

0 

■P2 

-Pj^/DY 

GI(I) 

URJ(I-l)  A2  UIJ(I-1)/2  > Al  VRJ(I-l)  A3  FRJ(I-l)/2  • 

I 

UIJ(I-l)  - A2  URJ(I-l)/2  -y  Al  VIJ(I-l)  -t  A3  FIJ(I-1)/2 
m3[ui2R(l,I)  ^ U12R(1,I-i3  -►  in4  [vl2R(l,  I)  Vl2R(l,I-l3- 

m3[ul2I(l,I)  ♦ Ul2I(l,I-l)j  m4  [vl2I(l,I)  V12I(1,  I-dJ - 

m4[ul2R(l,I)  -t  U12R(1,I-1)]  n4  |yl2R(l,I)  V12R(1 , 1-10 - 

m4  [U12I(1,I)  + U12I (1,1-1)]  4-  n4  [vl2I(l,I)  V12I(1, 1-l)]- 

(P3  - m2)[u12R(1,I)  4 U12R(1,I-1)]  + p4[vl2R(l,I)  ^ 

(P3  - M2)[ai2I(l,I)  •*.  U12I (1,1-1)]  ♦ P4[v12I(1,I) 
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A. 5 ( Continued ) ; 


3,  Matrix  equation  (riqht  hand  side  continued): 


- URJ(I-l)  - GR(I-1)/DY  - m2  GI(I-l) 

- IIIJ(I-I)  - GI(I-1)/DY  -►m2  GR(I-l) 

- VRJ(I-l)  - nj^  GR(I-1)/dy  - n2  GI(I-l) 

- VIJ(I-l)  - nj^  GI(I-1)/dy  -►  n2  GR(I-l) 

4-  V12R(I-1)]  4-  kVM2[F2I(I)  + F2I(I-1  ^ 4-^  M^URJCl-l  ) 
4*  V12I(I-i3  - k«M2[F2R(I)  4-  F2R(I-L)]4-rf  M^UIJd-l) 


-IJM^k  FIJCI-l) 
4-»M^k  FRJ(I-l) 


T>i  GR(I-1)/dY  - P2  GI(I-l) 
GI(I-1)/DY  4-  P2  GR(I-l) 


Blade-Two  Sxirface  Points. 


1.  References:  section  3.9.1;  figtire  3.9.1  A. 

2.  Auxiliary  constants:  X = X2(I,I)  “ X2(l,l) 

3.  Equations: 

V12R(I,I)  = - COS(^)  - k(X-b)SIN(^) 
V12I(I,I)  = - SIN(^)  - k(X-b)COS(yu() 
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A. 7 Equations  for  the  ui  Velocity  at  Blade -TVo  Surface 


Points. 

1.  References:  section  3,9.1;  figure  3,9.1  A. 

2.  Auxiliary  constants  (also  used  in  A.9  equations); 

DXA  = X2(I,J)  - X2(I-1,J)  XLA  = ^ ?i(I-l,j3 

DXB  = X2(I,J)  - X2(I,J-1)  XLB  = i[V(I,J)  4*  MI.J-D] 

VAR  = i|yl2R(I,J)  > V12R(I-1,J^  DLB  = >^I,J)  - X(I,J-1) 

VAI  = i[vi2I(I,J)  V12I(1,J-13  Al  = 2 k m2 

DVAR  = V12R(I,J)  - V12R(I-1,J)  A2  = k Al/2 

DVAI  t V12I(I,J)  - V12I(I-1,J) 

SLA  = i 7v(i-i,J)^tI  SLB  = X(i,J-l)*iI 

DLXB  = |>(I,J+1)  - X(I,J)]/[^(I,J-*.l)  - X2(I,j3 

3.  Matrix  equation: 

■I  ^ tiLXB  DXA  _ A2  DXA^  , Al  DXA 
4 XLA  • 2 XLA  2 XLA 

Al  DXA  , DLXB  DXA  A2  DXA^ 

2 XLA  ^ 4 XLA  “ 2 XLA 

U12R(I-l,J)[l  - U121(I-1,J) 

U121(I-1,J)  jl  - - U12R(I-1,  J)  ^ 


ljl2R(I,J) 

U12I(I,J) 


DVAR 

•►•STS' 


. A2  VAR  DXa2 
*••  XLA  SLA  •► 


F2R(I-1,J) 


DVAI  A2  VAI  DXa2 

XIA  §LA  •► 


F2I(I-1,J) 


A2  DXA 
XLA 

A2  DXA 
XLA 


4,  Ncte:  at  surface  points  I s J. 
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A«8  Eqiaations  for  at  Blade-Two  Surface  points. 

1.  References:  section  3,9.1;  figure  3,9,1  A. 

2.  AtDciliary  constants:  as  defined  in  A.7. 

3.  Equations: 

F2R(I,J)  = F2R(I-1,J)  4*  {i[ul2R(I,J)  4-  Ul2R(I-l,J)] 

4*  VAR/SLA^  DXA 

F2I(I,J)  = F2I(I-1,J)  4.  {iIul2I(I,J)  4 Ul2I(I-l,J)] 

4*  VAI/SLa}  DXA 

4.  Note;  at  surface  points  I s J, 

A, 9 Equations  for  the  Velocities  ui  and  at  Curved-Mesh 
Grid  Points. 

1.  References;  section  3.9.2;  figure  3.9,1  B. 

2.  Auxiliary  constants;  as  defined  in  A,7, 


3,  Matrix  equation: 


1 ^DLXB  DXA_A2  DXA^  _Al  DXA 

4 XLA  2 XLA  2 XLA 


Al  DXA  1 . . DLXB  DXA  _ A2  DXA^ 

2 XLA  ~ 4 XLA  ""  2 XLA 


1 A2_rora£  ^A1  DXB  1 
^ XLB  2 XLB  2 XLB  SLB 


Al  DXB  , , DLB  _ , A2  DXB^ 
2 XLB  4 XLB  2 XLB 


. 1 ^ A2  DXA^ 
SLA  i kLA  SLA 


1 — A2  DXa2 
SLA  2 XLA  SLA 


■X.  A2  DXB^ 

2 XLB  SLB 


_1 . A2  DXB^ 

SLB  2 XLB  SLB 


(Note:  continued  on  next  page.)  j 

i 
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A«9  (Continued): 

3,  Matrix  equation  (unknown  vector  and  right -hand- side 
matrix  continued): 


ui2i(i,j)  ui2i(i-i,j)  [j.  - ^ • 


U12R(I,J) 

U12I(I,J) 

V12R(I,J) 

V12I(I,J) 

[l  - 
[i  - 


A2 


4 XLB 


DLB  . ^ 
4 XLB  2 


DXb21  _ 

ias”j 


■V  U12I(I-1,J)  A3-.  PXA  V12R(I-1,J)  H-  -t-  A2  DXA^  "|  ^ 

’ 2 XLA  L SLA  2 XLA  SLaJ 

- U12R(I-1,  J)  V12I(I-1,J)  r,  i X A2  DXa2  ~|  ^ 

2 XLA  L SLA  2 XLA  SLaJ 

H>UI2I(I,J-1)  Al  DXB  ^V12R(I,J-1)  - A2  DXB^  ~|  ^ 

2 XT^  L SLB  2 XLB  SLbJ 

- U12R(I,  J-1)«A1_2^ -e  V12I(I,  J-1)  Hi-  ~ A2  DXb2  "I  ^ 

2 XLB  * L SLB  F5rf.B  5lbJ 


-►  F2R(I-1,J)  A2_J3^ 
XLA 


^ F2I(I-1,J)  A2  DXA 
XLA 


F2R(I, J-1)  A2  DXB 
XLB 


-I-  F2I(I,J-l)_A2_DJffi 
XLB 


A.IO  Equations  for  at  Curved-Mesh  Grid  Points. 

1.  References:  section  3.9,2;  figure  5.9.1  B. 

2.  Auxiliary  constants:  as  defined  in  A. 7. 

3.  Equations: 


P2R(I,J)  = ilF2R(I-l,J)  ^ F2R(I,J-1)] 

•V  |•^U12R(I,J)  4-  U12R(I-1,J)  ^ V12R(I,J)  V12R(I-1.  J) 

SLA  ■ 

•4*  i{ui2R(I,J)  -V  Ul2R(I,J-l)  - V12R(I.J)  - Vl2R( I.  J-l )‘ 

SLB  < 


DXA 

DXB 


F2I(I,J)  = i[F2I(I-l,J)  +F2I(I,J-1)] 

-►  i[ul2I(I,J)  ^ U12I(I-1,J)  Vl2I(I.J)j^Vl2I(I-l^J)^nYA 
4-  •5-^U12I(I,  J)  ^ U12I(I,  J-1)  - V12I(I.  J)  Vl2I(I.J-l)|nyp 
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CC^PUTi.NG  THE  CURVED  MESH  POINTS  (1,2),  AND  (2,2) 
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